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FORWARD 


This report describes work performed by Software and Engineer- 
ing Associates, Inc. (SEA) under NASA contract NAS8-35046, 
titled Improved Two-Dimensional Kinetics (TDK) Computer Program. 
The TDK computer program is a primary tool in applying the 
JANNAF liquid rocket engine performance prediction procedures. 
The purpose of this contract has been to improve the TDK com- 
puter program so that it can be applied to rocket engine designs 
of advanced type. In particular, future orbit transfer vehicles 
(OTV) will require rocket engines that operate at high expansion 
ratio, i.e. in excess of 200:1. Improvements to the TDK com- 
puter program are necessary if it is to be applied engines of 
the OTV type. This study represents a first step in extending 
the TDK computer program so that it can be used for these 
operating regimes. 

The method of approach used in this study has been to 
transfer technology from the Solid Performance Program (SPP) to 
the TDK computer program. The SPP was developed by SEA for the 
Air Force under contract F0il6n -77-C-0046 . In this manner, the 
TDK computer program has been equipped with a boundary layer 
analysis, and a limited capability for treating two-phase flow. 
In addition, a method has been added that will automatically 
calculate the effects of regnerative cooling on thrust chamber 
performance. 

Project manager for this project has been Mr. Gary R. Nic- 
kersor. Dr. Lanh Dang has performed the computer programming 
and has provided engineering support. The project was very much 
aided by the helpful support of the contract monitor, Mr. Klaus 
W. Gross, and by Mr. A. Krebsbach. 
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ABBREVIATIONS 


BLIMP Boundary Layer Integral Matrix Procedure, JANNAF 

ooiuputer pi‘Oj^raui» 

BLM Boundary Layer Module, computer program 

JANNAF Joint Army-Navy-NASA-Air Force 

M0C Method of Characteristics, module of SPP and TDK 

0DE One-Dimensional Equilibrium, moaule of SPP and TDK 

ODK One-Dimensional Kinetics, module of SPP and TDK 

PSS Performance Standardization Subcommittee of JANNAF 

SEA Software and Engineering Associates, Inc- 

SPP Solid Performance Program, JANNAF computer program 

TDE Two-Dimensional Equilibrium, module of TDK 

TDK Two-DlTensional Kinetics, JANNAF computer program 



1 . INTRODUCTION 


The Two-Dimensional Kinetics (TDK)^ computer program is a 

primary tool, in applying the JANNAF liquid rocket thrust chamber 

2 

performance prediction methodology. This computer program and 
the performance prediction methodology were originally developed 
under the auspices of the Performance Standardization Subcommit- 
tee (PSS) of the JANNAF. The goa . of the PSS is the development 
of a methodology that Includes all ■'spects of rocket engine per- 
formance from analytical calculation co test measurements, that 
is physically accurate and consistent, and that serves as an 
industry and government reference. 

Recent interest in rocket engines that operate at high 
expansion ratio, such as most OTV designs, has required an 
extension of the analytical methods used by the TDK computer 
program. The study reported here is a first attempt in 
extending TDK so that it is applicable to the operating 
conditions imposed by these high performance OTV designs. 

High performance can be obtained with large area ratio 
thrust chambers using oxygen with either hydrogen or hydrocarbon 
fuels. Due to the large area ratios, the combustion products 
are expanded to low pressures such that condensed phases of wa- 
ter may form in addition to soot particles with hydrocarbon fuel 
application. The existing TDK analytical nozzle flow model for 
liquid rocket engines did not possess the capability of particle 
flow simulation. However, a recently developed Solid motor 


Ref. 1. Nickerson, G. R., Coats, D. E., and Bartz, J. L. 

"The Two-Dimensional Kinetic (TDK) Reference Computer 
Program," Engineering and Programming Manual, Ultra- 
systems, Inc., December 1973» prepared for Contract No. 
NAS9-12652, NASA JSC. 

Ref. 2. "JANNAF Rocket Engine Performance Prediction and Calcu- 
lation Manual", CPIA 246, April 1975. 
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Performance Program (SPP) does contain the simulation of parti- 
cle flow . This capability, which assumes gas-particle 
equilibrium, has been transferred from the SPP to the TDK 
computer program. 

The engine specific impulse value includes boundary layer 
losses. At present, this loss is identified by feeding results 
from the TDK program solution to a boundary layer code which 
calculates friction losses along the wall and heat energy 
transferred through the wall. In a regeneratively cooled 
engine, most of the heat energy transferred to the wall is 
returned to the main combustion chamber in the form of increased 
propellant enthalpy. In addition, the friction effect requires 
a displacement of the wall contour in order to conserve mass 
flow through the nozzle. Both phenomena necessitate the 
recalculation of the two-dimensional flow field at least one 
time to achieve a mass flow and energy balance. At present, 
this process is conducted by iterating using two computer 
programs, the two-dimensional kinetics (TDK) and the Boundary 

4 

Layer Integral Matrix Procedure (BLIMP) program . Both of these 
programs are quite complex, especially the boundary layer code, 
and they must be run in series. 

For study purposes, a simplified approach is very useful, 
especially when many computer runs are required with only one or 
a few input parameters to be changed at a time. The 


Ref. 3- Nickerson, G.4R., Coats, D. E., Hermsen, R. W. , ”A Com- 
puter Program for the Prediction of Solid Propellant 
Rocket Motor Performance", Volumes I, II, III, Interim 
Technical Report, for Air Force Rocket Propulsion Labo- 
ratory, AFRPL-TR-80-3^, April 1981. 

Ref. 4. Evans, R. , "Boundary Layer Integral Matrix Procedure", 
BLIMP-J User's Manual, Aerotherm Division/Acurex Corpo- 
ration, Prepared for NASA George C. Marshall Space 
Flight Center under Contract NAS8-30930, July 1975. 
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Oriel-Dimensional Equilibrium (ODE) module, the One-Dimensional 
Kinetics (ODK) module, and the Two-Dimensional Equilibrium (TDE) 
module are such simplified options within the TDK computer 
programs The boundary layer program (BLIMP) also contains 
various simplified options. However, the program has a large 
computer core demand, needs long execution times, and the input 
is rather difficult to prepare. 

Recently, a new Boundary Layer Module (BLM) was developed 
for solid motor nozzles. This program holds all the essential 
analytical and operating features to produce an efficient coup- 
ling with the TDK program. As part of the study reported here, 
the BLM was taken from the SPF and coupled to the TDK program. 
The capability of treating regeneratively cooled nozzle walls 
was aaded to both BLM and TDK. 

The work carried out in this study is summarized below. 

a) The current method of condensed phase flow treatment in 
the ODE and ODK modules of the SPP program have been 
incorporated in the TDK program with appropriate 
changes. Specific attention has been given to the 
presence of water droplets and ice crystals and of 
solid carbon particulates. 

b) Using subroutines existing in the SPP, the ODE module 
of TDK has been modified so that it will compute trans- 
port properties for the nozzle exhaust composition. 
This information is required by the BLM. Numerous other 
features contained in the SPP, such as wall contour 
options, have also been transferred to TDK. 

c) The ODE module of TDK has been modified so that it will 
provide tables of gas properties versus temperature for 
the stream tube mixture ratio adjacent to the nozzle 
wall. This information is required by the BLM. 
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d) The Boundary Layer Module (BLM; of the SPP has been 
coupled with the TDK program such that It will auto- 
fuatically calculate the wall boundary layer after 
execution of the TDS or TDK modules, if so desired. The 
resulting wall displacement and the heat transferred 
tnrough the wall increasing the energy level of the 
cooling propellant are used to repeat a TDE or TDK 
calculation once with the appropriate parameters 
(wall contour and propellant enthalpy) updated. 

e) The boundary layer module has been equipped with a 
feature to account for the heat transfer along wall 
contour sections where regenerate cooling takes place. 
The wall heat flux, as determined by the BLM, is 
integrated and the result (multipled by an input 
efficiency factor) is returned to the main combustion 
chamber in the form of increased propellant enthalpy. 

f) Two proof cases have been run on the revised program 
for the purpose of verification. The first case is a 
Pratt 4 Whitney RL-10 engine configured to operate at 
an expansion ratio of 205. The second case is a 
Rocketdyne Advanced Space Engine (ASE) configured to 
operate at an expansion ratio of 400. Performance 
measurements exist for both of these engines. 

g) A User’s Manual has been prepared for the improved TDK 
progi’am. This has been done by revising Section 6 of 
Reference 1. A copy of this User’s Manual is attached 
as Appendix A. 

h) The Improved TDK computer program has been delivered to 
NASA-MSFC for use on the UNIVAC 1100/82 computer. The 
program is also operational at SEA on the DEC VAX/780 
computer . 

Technical details of the above work items are reported in 
Sections 2 and 3- 
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2. METHOD OF APPROACH. 


Prior to the work performed durinj? this contract . two 
versions of the TDK program have been generally available. The 
first of these is the version documented in Reference 1 . A 
number of improvements have been made to this program by SEA, 
such as those described in Reference 5. However, the basic 
program is as presented in Reference 1 . A second version of the 
TDK program has been developed by NASA-MSFC and has been 
distributed by them. This version of TDK contains many 
programming modifications. For the most part, these 
modifications are improvements, to the internal structure of the 
program and its input and output features. The basic analytical 
models have not been changed. The program is undocumented. 

For this contract it was necessary to combine the beat 
features of both of tne above programs, fo’^ming a "baseline" 
computer program, and to then transfer advanced features 
contained in SPP to this baseline program. Most of the features 
of interest are contained in the ODE and ODK modules of SPP and 
concern fluid chemistry. 

Features of SPP other than the fluid chemistry were also 
transferred to TDK. In particular, the wall geometry options in 
SPP were installed into TDK. First, all of the wall geometry 
options were moved forward in the code and unified in one 


Ref. 5. Nickerson, 0. R., "Modification to the Two-Di...ensional 
Kinetics (TDK) Computer Program", Software and Eng- 
ineering Associates, Inc., prepared for NASA George C. 
Marshall Space Flight Center, under Contract No. 
NAS8-31988, October 1976. 
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Namelist input data set, called $DATA. Next, the spllhe fit 
option for ODK was modified so that the DDK printout could he 
requested at input values of area ratio. The option for a cone 
with specified end point was added to ODK. An option for a wall 
defined by a skewed-parabola was transferred from SPP and added 
to both ODK and TDK. These geometry options are described in 
the revised User's Manual, Section 6.3.3 (see Appendix A). 

The "problem card" feature of TDK was replaced with « 
"module selection" feature of the type used by SPP, i.e., module 
selection by using input flags. This method a] lows greater 
flexibility than the "problem card" method, and is described in 
the revised User's Manual, Sec'- ion 6.3.1. 

Wi :h the completion of the transfer of technology from SPP, 
the condensed flow capability and the transport property 
calculations items l.a and l.b) were essentially complete. A 
discussion of the gas-particle properties now used in TDK is 
presented in Section 3.1. 

The requirements for coupling BLM with TDK go well beyond 
the usage of BLM in SPP. SPP does not require a second pass 
through the Method of Characteristics (MOC) solution, and there 
is no regen cooling to consider. These two topics, BLM and 
regen cooling, are discussed in Section 3.2 and 3.3. 

Further work is being carried out on TDK by SEA under 
contract NAS8-3i*974. The work being done is primarily concerned 
with adding a shock-wave capability to TDK. Work is being 
coordinated between the two contracts so that one final computer 
program will result. 
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3. TECHNICAL DISCUSSION. 

The technical effort i3 discussed under the three topics 
given he low. 

Fluid Properties, 

The Boundary Layer Module, and 
Regenerative Cooling. 


3.1 Fluid Properties . 

In order to expand the fluid property capabilities of TDK, 
changes of three basic types were made to the program, as 
follows : 


a) The calculation of transport properties (gas viscosity, 
conductivity, and Prandtl number), 

b) The calculation of gas properties adjacent to the 
nozzle wall, and 

c) Core flow gas-particle equilibrium chemistry for the 
ODK, TDE, and TDK modules. 

Each of thes' items is discussed in the following subsections. 
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3.1a Near Wall Gas Transport Properties. 

Values of gas vissosity, y , gas conductivity, < , PrandtJ 
number, ?r, are required by BLM. These transport properties arc 
evaluated along the wall streamline as a function of temperature 
using the ODE module. It is assummed that the expansion is in 
equilibrium, including condensed species. In the model these 
properties are then computed as "frozen” properties using the 
local gas composition with condensed phases included. 


The viscosity and thermal conductivity of the individual gaseous 
species are calculated from formulas given in Reference 6 as 
follows : 




4.15822 X 10 




ORIGINAL PAGE (S 
OF POOR QUALITY 




“^i ‘ M, 


w , 


■<•‘'5 * '-32 


W . 
1 


The required Leonard Jones parameters, are internally 

stored in the computer program for 206 gas phase species. (See 
Subroutine MUK, Volume II, of Reference 3, for a list of the 
species.) The viscosity of the mixture is calculated from 

7 

Wilke’s semi- mpirical formula . 





(1 


N 


^ij X. ^ 



where N is the number of species, x^ the mole fraction of 
species i, and 0^^ is defined by: 


Ref. 6. Svehla, R. A., "Estimated Viscosities and Thermal Con- 
ductivities of Gases at High Temperatures," NASA 
TR-132, 1962. 

Ref. 7. Bird, R. B., Stewart, U. E., and Lightfoot, E. N., 
Transport Phenomena , John Wiley i Sons, i960. 
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The results calculated as described above are used to 
provide the BLM with properties usinf the following procedure. 


The viscosity at the nozzle throat, u*, and the 
corresponding temperature, T*, are calculated. Using these 
values, an exponant,<«) , is computed such that the expression 

M = y* (T/T*)“ 

provides a best fit in a least square manner to the chamber and 
e.:it values for viscosity. Next, a constant value for Prandtl 
number is computed as 



where y , k, , and C„ are evaluated at the nozzle throat. The 

P 

"frozen” Cp value is used. The values transmitted to BLM are 

then u*, T*, to , and Pr^*. Actually, BLM only uses the ratio 

Pr/y, and does not require k or the frozen versus T. 

P 


Ref. 8. Mason, E. A., and Saxena, S. C., Physics of Fluids, 
Vol. 1, No. 5, pp. 361-369, 1958. 


a 



3.1b 


Near Wall Gas-Particle Mixture Properties . 


Values of specific heat at constant pressure, Cp, and ratio 
of specific heats, Y , static enthalpy, n , and the gas constant 
R, are required by BLM. These properties are evaluated along 
the wall streamline as a function of temperature using the ODE 
module. The tables are prepared using a series of (T,S) 
equilibrium calculations, where T varies from 600°R to 7000°R at 
200'^R increments. The chamber entropy value is used for S. 
Values at 100°R are then extrapolated and added to the table. 
The table is printed with the BLM output. An example of this 
output is presented in Figure 1 . 

Since the ODE module did not contain a (T,S) option, it was 

adcied using the following procedure. Using known values of 

T,, and Y_, and the given value of temperature, T, a first 
a a ^ 

estimate for pressure P was found as 
P^ ^ = P^ (T/T,). ^ ® 

a a 

The (T,P) option of ODE is used to obtain as a function of 

p(i)^ The procedure is iterated using the secant method 
(subroutine ITER plus a driver) to find P^^^ such that 

|(S^^^-S^)/Sq |< 5 E-5 

where S is the chamber entropy. The procedure is repeated for 
c 

each temperature in these tables. The procedure is internal to 
the program and not callable through the ODE input. 

If the chemical system contains a condensed phase, then 
equilibrium solutions made during phase change can yield values 
for Y that are unacceptable to BLM, e.g., values of Y < 1. 


10 



GAS PRUPe’^TIES tables 




OR/OTML F;»Ge rs 

OF POOR QUALITY 


s: 

02 


9990900099 ^ 

99000909904 

«i»l 

L»fV^<V 4 V<^ 4 rwiM<v<Mnt< 

« •••••••••• 


ni nii^ 

A1 M A| 

Al 

ni Ai Ai 

Ai 

“M 

Ai 

A# 

AI 

AI 

Ai 

At AI 

A 

000 

0 0 9 

0 

e 0 0 

0 

0 

0 

0 

0 

0 

0 

0 0 

0 

♦ ♦ 

♦ ♦ ♦ 

♦> 

♦ ♦ ♦ 

♦ 

♦ 

4> 

♦ 

♦ 

♦ 

♦ 

♦ ♦ 

♦ 

19 19 Mi 

UJ 19 19 

19 

19 19 19 

W9 

19 

Mi 

19 

19 

19 

19 

Mi 19 

19 

0 0 

Ai 9 

0 

9 ^ rg 

9 

9 

0 

0 

9 

Ai 

9 

9 n 

A* 

^ 9 «n 

9 9 C 

9 

A| 9 9 

in 

Ai 

•m 

0 

« 

9 

0 

9 in 

9 

•*►0 9 

«u9l 0 

mm 

^ 9 nj 

Ai 

•A 

9 

0 

« 

0 

AI 

iA 9 

9 

IV Ai rw 

»n in 9» 

9 

9 in 0 


« 

9 

0 


fA 

9 

in 0 


• * • 

« • • 

• 

# • • 

• 

• 

• 

• 

• 

• 

• 

• • 

« 

«««»«■ 



— 

— 

— 

•* 

rj 

Ai 

AI Ai 

Ai Ai 

Ai 


o; 

4^ 


<D 

n 


c 

C-, 


CO 

:3 

CO 

(V 

> 

QC 

'O 

5 


oooooooooooooo oooooooooooooooooooo 

tei Mi Ml «kl (M lu Ml( 

X o o < 



••••••••••••♦• 

90VV99999999VQ 99M 09999009990999999 
00000000000090 90000000000000000000 





|A|IWfW«MA|IMP««l«fi 4 % 999 l/>l/>tf >0000 


oooooooooooooo 


oooooooooooooooooooo 

oooooooooooooooooooo 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 


^9 ooo^oini^i^iKo (Vlur^ooli^o•««n^»#vfMlWoo^N>olhl^i•« 

ino«Do 9 i««P*»M»n 9 o#*> — 

otnov^ooomo 9 inr^r^i*^ >»» 9 »» 09 o^»«ir»in«oooo<MAi«»f'-© 
«Hr^909»^^O9«O^I9 4^9OP^O*«9P«O9Ol«l«0K»999 9«P<*I 

tfioii«^Hi 9 yp>of^oo 9 «w«m«T 9 tf%mo ok#^o 


o 

#1 

2 si 


o o in 

0 o o 

1 I » 
^ Ul 4 l 

a m 
3 •- Oi o 

X O O 9 
09 •• 


m m 
o o 

• • 


o 

o 


0 o < 

1 • 

Ml -Jj « 

<M P 

•- O I 
M Kl 

«o o 


I m in I 
10 0 4 

• f 


P* 9 
9 m 
fM 9 


u> u> 
o o 




in tn 
o o 

t t 

iAi UJ 

9 in 

9 

9 « 
O O O «• «»| 


in 


9^ m 

0 o 

1 f 

nAt <AI 
ni 9 

o »n 
fw — 
o o 


m in 

0 o 

1 I 

UJ UJ 
fst 9 

o 

9 9 

in 


in tn 
o o 
t t 
Ml UJ 
9 9 

in *n 
o c* 

«» m 


in in m in m m m 

0 o o o o o o 

1 t » I t I f 

Ul tel uJ Mi Mi Ji Mi 
0 0 9 O ^ 4 J 

9 9 c nj 9 in tn 
tn *n c o m M 9 

9 O O 9 *n 9 


m in 
e o 
I • 

UJ Mi 
O « 

m *n 
o m 
o o 


mo«*«*w*««mnininim»npnpn ^>999999 otntnininmmooooo^ 


0 

0 

0 

0 

0 

c 

0 

« 

0 

0 

0 

fM 

A- 

»u 





A> 

fM 


A» 

A«. 

A> 

p«*. 


A» 

A>. 

0 

0 

0m> 

A* 

A* 



0 

0 

0 

C 

0 

0 

0 

C 

0 

c 

0 

0 

0 

0 

0 

0 

0 

c 

0 

0 

0 

© 

0 

0 

0 

0 

0 0 

0 

0 

C 

0 

0 

♦ 

M 

♦ 

♦ 

♦ 

•♦■ 

U' 

♦ 

♦ 


♦ 

♦ 


♦ 

M 

♦■ 


♦ 

♦ 


♦ 

♦ 

J> 

•A 

♦ 

♦ 

♦ 

M 


♦ 


A 

-A 

A 

Ml 

19 

19 

Ui 

19 

19 

19 

Ui 

19 

UJ 

19 

Mi 

Ui 

Ui 

Ui 

19 

UJ 

19 

Ui 

Ui 

Ui 

19 

19 

19 

19 

19 

19 

Ui 

Ui 

19 

Ui 

Ui 

Ui 

19 

0 

Ai 

9 

9 

A* 

in 

9 

lA 

AI 

e 

9 

« 

X 

in 

0 

9 

AI 

0 

e 


A* 

m 

0 

e 

9 

Ai 

^9 

9 

#0 

« 

0 

0 

91 

AI 

tA f>» 

£ 

0 

•A 

0 

0 



0 

0 

0 

IM 

fu 

0 

9 

0 

0 

0 

9 

9 

< 

AI 

»A 

9 

© 

X 

n 

tn in 

9 

A* 

A- 

a 

A. 

X Ai 

9 


0 

n 


C 

9 

0 

n.> 

tn 

e 

9 

r«» 

9 

•• 

in 

9 

0 

9 

in 

9 

9 

ii> 

9 

m 

in 

9 

A# 

0 

0 

tn 

A* 

■* 

9 

4A 

tn 

AI 

Ai 

AI 




0 

© 

9 

9 

c 

0 

Ai 

r- 

AI 

c 

Ai 

Ak 

M 

«n 

9 

Ai 

in 

0 

0 


9 

•A 


9 

0 

• 

• 

0 

• 


• 


• 

• 

• 

• 

9 9 

m 

9 

to 

C 


A* 

4)oinin94A9iAt 


— 

»A 

9 


aI 

Ai 

"A 


( 1 ) 


AI 

Ai 

AI tJ> 

*A 


m 

tA 

iA 

♦A tA 

*A 

•A 

•A 

•A 

•A 

iA 

lA 

tA 

•A 

»A 

tA 

«A 

tA 

lA 

•A 

«A 

•A 

iA 

iA 


•A 

PA 

•A 

e 

C 

C 

© 

0 

0 

0 

0 

9 

0 0 

0 

0 

0 

0 

0 

0 

O 

e 

0 

0 

9 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A 

♦ 

♦ 

M 

♦ 

A 

A 

A 

A 

A A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

A 

Ui 

19 

Mi Ui 

19 

19 

19 

Ui 

Ui 

UJ Mi 

19 

19 

19 

Ui 

Mi 

19 

Mi 

Ui 

UJ 

19 

19 

Ui 

19 

19 

19 

Ui 

19 

Mi 

Ui 

19 

Ui 

UJ 

19 

0 

0 

0 

0 

0 

0 

e 

0 

0 

0 0 

e 

0 

0 

0 

0 

0 

0 

0 

0 

e 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

C 

0 

e 

0 

0 

0 

0 

e 

0 

0 

c 0 

0 

e 

0 

e 

e 

e 

0 

e 

0 

e 

e 

0 

0 

0 

0 

e 

0 

0 

c 

0 

0 

0 

© 


0 

0 

0 

0 

0 

0 

e 

0 

0 0 

0 

© 

0 

e 

0 

0 

0 

0 

0 

e 

0 

0 

0 

C 

w 

0 

9 

0 

© 

9 

0 

0 

© 

9 

0 

0 

0 

Ai 

T 

0 

c 

0 

AI a 

0 

c 

e 

Ai 

a 


« 

0 

Ai 

a 

0 

« 

0 

Ai 

a 

0 

C 

9 

AI 

a 

0 

« 

c 

• 

^ m 


\S 

• 

' 0 
«» 

• 

A, AI Ai' Ai Ai' iA 

lA A> tA fA a 

a 

a 

a 

a in m tn in in 0 

0 0 

0 

0 



t- 


•*»v:pn 9 vrOP^«. 9 ©«-ni »^9 


lno ^»99 o*»»m '^9 no ^®9 o^-n'^na 
«.«.««^«Miv'vi'MfV'vrvntmnirMfnpn 4 npn*n 


u 


i) 



Since condensation is an isothermal process, using (T,S) \ 

solutions avoid this problem. It also avoids the "temperature 

out of range" problem that can occur had other, existing options 1 

K A ^ .-HI 

uociu* 


The fluid properties shown in ""able 1 are h, u , < , Cp,Y , 
and R versus T°R. The values of u and < are prepared as 
described in Section The values for h, Cp, y , and R are 

prepared as described above, i.e., they are equilibrium 
properties Including condensed phases (no gas particle lag). 



3.1c 


Core Flow Gas-Particle Equilibrium Chemistry. 


The TDK program described in Reference 1 contains only the 
fliosL rudimentary capability to compute the effects of two-phase 
flow (gas plus condensed phases) on nozzle performance. In that 
program only the ODE module is capable of considering the 
existence of different phases. This is done by treating the 
entire chemical system as if it were in chemical equilibrium and 
the particulate and gas phases as if they were in dynamic and 
thermal equilibrium. For high pressure systems and where the 
characteristic particle sizes are small, the above assumptions 
are quite reasonable. For systems where there is a significant 
departure from either of the above assumptions, then the TDK 
code described in Reference 1 does not have the capability to 
even assess the importance of these effects, let alone to 
quantitatively predict the magnitude of them. Systems where 
two-phase flow phenomena could be of importance are those with 
soot formation using hydrocarbon propellants, and with the 
formation of water droplets and ice in low pressure, high 
expansion ratio H 2 /O 2 systems. 

In order to improve the capability of TDK with respect to 
two-phase flow, the ODK module of SPP was installed in TDK. The 
chemistry subroutines of ODK also serve the same function for 
TDK, so that TDK has been given the same two-phase flow (no lag) 
capability as ODK. Unlike TDK, the SPP has no two-dimensional 
kinetics capability. 

The Two-Dimensional Equilibrium (TDE) option has been 
updated to consider two-phase flow in the same manner as the ODE 
module, but for a two-dimensional calculation. Source flow 
calculations were carried out to verify the program. Under this 
assumption TDE was shown to calculate the same results as ODE 
for all streamlines. The TDE capability allows the user to 
assess the effect of particles on specific impulse, and also to 
indicate the possibility of the onset of condensation. Further 
development alor.g these lines would require the generation of 
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kinetic models to predict both condensation and mass transfer 
between phases and is beyond the scope of this study. Implicit 
in the above is the assumption that gas-particle lag effects are 
small. If this assumption is violated then the SPP code should 
be used to estimate the performance loss because the two-phase 
flow loss becomes the major thrust loss in the system. 



3.2 


The Boundary Layer Module (3LM) 


The BLM was developed for SEA by CBC Enterprises. 

Documentation for the BLM as delivered to SEA is presented in 
Reference 9. The purpose of the BLM is to compute compressible 
laminar and turbulent flows in axisymmetric nozzles. It can 
also be used to calculate two-dimensional and axisymmetric 
external flows. The method utilizes an efficient two-point 

finite difference method developed by Keller and Cebeci^^. 
Turbulence modeling is achieved through the use of the 
Cebeci-Smith eddy-viscosity formulation^^ which has been tested 
for a large class of flows with various boundary conditions. 

The BLM has been modified by SEA to interface automatically 
with the SPP and with the TDK programs. TDK supplies all of 
the inputs required by BLM. 

In order to be able to use BLM for large expansion ratio 
nozzles, the boundary layer calculation was divided into 
adjacent axial segments, and the calculation restarted 

automatically at each segment. This feature is also useful when 
nozzles with regen cooling are being analyzed. A segment can be 
placed at the start (or end) of the cooling circuit where there 
is a discontinuity in the wall temperature (or heat flux), and 
the calculations can be restarted with an integration step size 
change. 


Ref. 9. Cebecl, T., "Boundary Layer Analysis Module", prepared 
for Software Engineering Associates, Inc., January 
1982. 

Ref. 10. Keller, H. B., and Cebeci, T. , "Accurate Numerical Met- 
hods for Boundary Layer Flows, Part 2, Two-Dimensional 
Turbulent Flows", AIAA J., 10, 1972, p. 1193. 

Ref. 11. Cebeci, T., and Smith, A. M. 0., Analysis of Turbulent 
Boundary Layers, Academic Press, N. ^. , 1 974 . 



Ac high expansion ratio , it is possible for the gas density 
in the boundary layer to become very low. If the mean free path 
in the boundary layer approximates the wall roughness height, 
the non-slip condition at the v/all can break down. To aid In 
the assessment of this condition, the mean free path, \ , at the 
edge of the boundary layer is calculated and printed, where 

A = 1 .26 U/ ( P\/rT ) (ft) 

At the end of the nozzle, the Knudsen number, Kn, based on 
displacement thickness, §*, is printed so that the approach to 
rarified flow can be monitored, where 

Kn = ^ / s! 

e 

The following options for plotted output have been added to 
the BLM by SEA: 

1) Momentum thickness, 6, versus axial position. 

2) Displacement thickness, <s*, versus axial position. 

3) Wall temperature, T^, versus axial position 

4) Velocity profiles at specified are ratios, or at 
specified axial locations. 

5) Temperature profiles at specified area ratios, or at 
specified axial locations. 

The analysis used by the BLM is described in Reference 9. 
Usage of the program within TDK is described in Section 6.8 of 
the User's Manual (see Appendix A), as follows: 

6.8.1 BLM Gas Properties 

6.8.2 Boundary Layer Edge Conditions 

6.8.3 Integration Step Size Control 

6.8.4 Regenerative Cooling Heat Transfer 

6.8.5 BLM Plotted Output 
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3.3 


Regenerative Cooling 


When an engine is cooled using a regenerative device, 

propellant (usually fuel) is routed around tne nozzle so that 

* 

heat is transmitted from the boundary layer to the coolant . 
This heat is then returned to the combustion chamber in the form 
of increased propellant enthalpy. When the nozzle wall 
temperatures are assumed known (IHFLAG = 0 option), the BLM will 
calculate the heat flux from the boundary layer to the wall, 

ft^(BTU/ft -sec). These values can also be input directly 
(IHFLAG = 1 option). If the coolant circuit extends from 

position Xq to position x^ and operates at an efficiency, n , 
then the propellant enthalpy entering the chamber will be 
increased by 

X 

e 

AH = “ n f n dA (BTU/lbm) 

W 

T *0 

Where 

is the total engine mass flow rate 

dA is the nozzle surface area differential, 2Trrdx 

Using the method outlined above, the BLM will compute increments 
of propellant enthalpy for up to 3 fuel or oxidizer circuits and 
print out the resultant enthalpy increments. These can be added 
to a later computer run by using the DELHI (1) input array. If 
BLM is to be automatically rerun, then the enthalpy increments 
will automatically be stored into DELHI (1) for the second pass 
through ODE, ODK, and TDK. These enthalpy increments can be 
calculated in two ways. If the enthalpy increase is distributed 
equally throughout the chamber, then 


Heat c.tn be extraotod from the boundary layer as a function 
of length (BLM analysis). The core flow (TDK analysis) does 
not include heat transfer effects other than to modify the 

f 

initial value for chamber enthalpy used by TDK . 
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AH, = AH 


i . e . 


DELHI (1 ) = AH, AH, etc. 

The second method is to assume that a fuel circuit adds 
enthalpy only to fuel, and an ox circuit adds enthalpy only to 
ox. It follows that 

for a fuel circuit 
r+1 


and for an ox circuit 

r.l ''l 

= <F7TT> - ‘W 

where r^^ is the mixture ratio of zone i and r is the overall 
chamber mixture ratio. 

For either method, the steady state engine cycle balance can 
be approximated as follows. First, calculate "adjusted tank 
enthalpies" for the fuel and for the oxidizer and input these on 
the reactant cards. These values must approximate the energy 
content of the propellant entering the main combustion chamber 
accounting for all energy gains and losses, except heat returned 
to the main combustion chamber by the regen cooling circuit(s). 
Estimates for these amounts ^re to be entered using the DELHI (1) 
input array. An estimate o* zero is usually satisfactory. Cor- 
rected estimates will be calculated by BLM and stored in 
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DELHI (1) for a second pass through TDK (or TDE). A second pass 
using these values will be executed automatically if IRPEAT = 1 
was input in the $DATA namelist. 

A second pass can be made without adjusting the enthalpy by 
setting IRPEAT = 2. In this case neither the ODE nor ODK runs 
will be repeated. The nozzle wall, however, will be displaced 
for TDK (or TDE). This option is useful if an adiabatic wall 
(zero heat flux) case is run. 

The boundary layer calculations are automatically started at 
the downstream end of the combustion chamber (position ECRAT of 
Figure 6-1, Appendix A). However, an option has been provided 
that allows the user to input values for the upstream boundary 
layer edge conditions and then start the boundary layer further 
upstream (see item lOFF, Section 6.3.2, Appendix A). This 
option can be important when calculating heat trrr^fer for the 
chamber-throat wall assembly. 
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4. 


CONCLUSIONS. 


Work items a) through h) described in Section 1 of this 
report have been corapletod: This has resulted in the teohnioal 

improvements to the TDK computer program that are discussed 
below. 

TDK has been equipped with a limited capability (zero 
velocity and thermal lag) for treating rocket nozzle exhaust 
flow containing condensed phase species. This has been 

•5 

accomplished by a transfer of technology from the SP?"^ to TDK. 

9 

A boundary layer program, known as BLM , has been installed 

in TDK as a module. The BLM was developed for SPP in response 

to requirements set forth by experts in the field at the 2nd SPP 
1 2 

Workshop . Further development of this module has been carried 
out by SEA to make it applicable to high area ratio liquid 
rocket engines with regen cooling. 

TDK has been modified to incorporate a capability for cal- 
culating the effect of regen cooling on engine performance. The 
nozzle wall temperature (or heat flux) profile must be specified 
by the user. The engine cycle balance must also be performed 
external to the program. 

The CPIA 246 standard procedure for computing engine 

2 9 

performance has been automated within TDK; but with BLM 

4 * 

substituted for BLIMP . A complete TDK-BLM-TDK calculation can 
be carried out in one computer run with a minimum of user input. 

Two sample cases were run for the purpose of demonstrating 
the program. These are 1) the Pratt and Whitney RL10 engine 
operating at an area ratio of 205.03:1 and 2) the Rocketdyne 




Advanced Space Engine (ASE), operating at an area ratio of 
400.71:1. Results of these computer runs have been sent to 
NASA-MSFC for evaluation. 

The improved TDK program has been delivered to NASA-MSFC 
converted for use on the Univac 1100/80 computer using the FTN 
compiler. The program is available from SEA for the DEC VAX 780 
computer c 

A User's Manual for these programs are presented in Appendix 

A. 
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5. 


RECOMMENDATIONS. 


Recommendations concerning the following topics are dis- 
oussed below: chemistry, low density, flow effects, test c?ses. 
input/output, and documentation. No recommendations are m..de 
here regarding shocK waves since TDK is currently undergoing 
modification for this purpose under Contract NAS£-3^974. 

m 

5.1 Chemistry. 


The TDK program contains two chemistry modules, which are 1) 
the chemical equilibrium solution subroutines of ODE, and 2) the 
streamline chemical kinetics subroutines of ODK. These latter 
subroutines are also used in a modified form to compute kinetics 
in a TDK calculation. No improvements to the CDE chemistry are 
recommended. However, several important improvements are 
recommended for the ODK (and TDK) chemistry. 

ODK should be modified to include input of "global" reac- 
tions and "hetrogeneous" reactions. For example, hydrocarbon 
decomposition can be characterized by reactions such as: 

^ CHm H 2 , other fragments. 

In order to effectively model the chemistry of propellants such 
as LOX/RP-1 , it is often necessary to include reactions of this 
type. Further examples are the conderisat<on of water, and also 
the kinetics of so^o formation. These processes can be 
Important in rocket engines. Further work is needed which will 
allow ODK to be used to assess the effects of these proceesses 
on rocket engine performance. 

A more advanced reaction rate processc-’ is available in the 
GKAP computer program. It is recommended that tnls processor be 
installed in ODK. This will add the following capabilities to 
ODK; 
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a) Improved treatment of third body type reactions (rate 
ratios assigned with each reaction). 

b) Input of reaction rates of the Landau-Teller form 
(necessary for r-)st T-V and V-V reactions). 

c) radiative reactions, such as 

HF(1 ) — HF(0) + PH0T0N 

The new rate processor will provide a firm basis for 
expanding the chemical kinetics capability of ODK. 

Other recommendations in the chemistry area are: 1) pro- 
viding thermodynamic data for the temperatures below 100°K, 2) 

replacing the thermodynamic data library with the latest data 
available from the NASA LRC, 3) updating of the kinetic rate 
data for reactions in the CFHNO systems that are important in 
rocket performance calculations. These items should be done on a 
"first cut" ’ '•sis and then reviewed by the JANNAF membership. 



5.2 


Low Density rlow Effects . 


Exhaust nozzle designs for future space engines may operate 
at expansion r-.tios so high that low density flow effects must 
be considered. If so, the following questions must be examined: 

a) At what conditions do low density effects at the nozzle 
wall require that a slip condition be imposed? 

b) Where does transition to free molecular flow cause a 
breakdown of the continuum flow assumption? 

c) Is the boundary layer interaction with the core flow 
sufficiently weak so that the JANNAF procedures can be 
used with confidence? 

Questions a) and b) above, can be answered in part by examin- 
ing the mean free path based on the average gas molecular 
weight. This item has been added to the BLM output along the 
boundary layer edge. For the boundary layer, question b) can be 
answered in part by examining the Knudsen number based on 
displacement thickness. This item has also been added to the 
BLM output. 

Question c) is more difficult in that it requires an 
analysis of the boundary layer-core flow interaction. It is 
recommended that an analysis of this interaction be conducted. 


2 ^ 



5.3 


Test Cases. 


The JANNAr standard performance prediction procedure was 
applied to two engines (the RL10 and the ASE) using the improved 

4 

TDK program. The BLM was used rather than BLIMP-J . Input data 
for the cases. Including data on the engine cycle balance, was 
provided by NASA/MSFC. 

It is recommended that further cases be run for the purpose 
of obtaining estimates of the accuracy of ODK, TDK, and BLM as a 
function of items that control the numerical methods. Cases 
should be run to determine if all of the program options work 
properly. 

Cases should also be run for the purpose of determining the 
role of the boundary layer in performance prediction. The 
following questions need to be addressed: 

1) How do the BLM predictions compare with BLIMP-J? 

2) How do BLM calculations with regen cooling compare with 
adiabatic wall calculations? 

3) What is the magnitude of performance gains due to 
"enthalpy pumping" in regen cooled engines? 
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Input and Output for 7DK . 


5% 4 


The following recommendations are made with respect to input 
and output for TDK and it-.? modules. 

a) A summary page should be added at the end of the output 
which lists in a clear format all of the key 
performance results. This page would correspond to the 
summary page listed by the SPP. 

b) Some, or all, of the following computer plots should be 
added to the program: 

Igp ^‘wall " 

®wall’ '"wall ^®^3us 

P,V,M, versus and 

P and M contour plots 

P,V,M, versus r at X .. 

6 X 1 0 

MOC plots for a window 

c) The linkage data files of ST P have been installed in 
TDK (see Reference 3» Volume III, Section 2.9). This 
capability allows a module of TDK to be run from a file 
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generated on a previous run. To be completely 
effective the file structure and programming must be 
modified specifically for TDK, and documented. This is 
recommended. 

d) The BLM output should be reviewed for the purpose of 
improving the format and to determine which quantities 
should be added, or deleted. 
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5.5 


Documentation. 


A complete set of manuals for TDK should be written 


discuss the following topics: 


1) Methods of analysis (all modules). 

2) Program structure and subroutine write-ups. 

3) User's Manual. 


which 


The documents should be written so that new users of the 
computer program can be systematically guided through the prin- 
ciples involved in the analysis. The analytical procedures 
should be related to the thrust chamber performance prediction 
methodology defined in CPIA 246, Reference 2. Test cases should 
be Included that demonstrate the differen.. f'^'^blem types solved 
by the program. Text shoulu be included describing each test 
case; its purpose, the engine type and c laracteristics, unique 
features, interpretation of output "esulcs, etc. 
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6.0 


PROGRAM USER'S MANUAL. 




Description of the Computer Progam Input. 


Ihe TDK computer program consists of five modules, ODE, ODK, TRANS, MOC, and 
BLM. All of these modules are required to perform a complete two dimensional 
non-equilibrium nozzle performance calculation with a boundary layer. Various 
options exist in the program, however, \diich exercise the above modules alone, in 
part, or in combination. 

Data is read by the program sequentially in the order required for the execu- 
tion of the modules. Ihls order is as follows: 


Thermodynamic data. 

Data common to the modules, 

(XE module inputs, 

ODK module ir^suts, 

IRAi^S module ir^uts, 

MOC module ir^uts, and 
BLM module iiq>uts. 

A more detailed description of these irput data sets is presented in Table b-1 . 
The documentation in which each of the data sets is completely described is also 
indicated in Table 6-1 . 

Of the data sets listed in Table 6-1 , only the $DATA data set is required for 
every computer run. Input of the other data sets is required only if the options 
they contain are to be used. For example, input of the thermodynamic data is not 
required if an existing thermodynamic data file is to be used. 


1 

( 

i 
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Table 6-1. Input Data Set Description 


Card Input 

Section No., 

Description 

THERMO 

« 

6.1 , 

thermodynamic data cards, see Tables 6-2 and 6-3 

• 

END 

UW T CPHS 

• 

6.1 .1, 

extension of thermodynamic data to temperatures 
below 300°K, see Tbble 6-5 

• 

END L0¥ T CPHS 

TITLE 

6.2, 

one or more title cards 

DATA 

6.3, 

data directive card 

$DATA 

SEND 

• 

VO 

data namelist for module selection, and ^cmetiry 

REACTANTS 

6.4, 

reactants directive card 

• 

6.4*1 , 

reactants cards, see Tables 2-6 and 2-7 

• 

6.4.1, 

blardc card required to end reactants cards 

OMIT 

6.4.2, 

cards to omit species 

INSERT 

6.4.2, 

cards to Insert species 

NAMELISTS 

6.4.3, 

0DE directive card 

$ODE 

6.4.3, 

0CE namelist 

SPECIES 

• 

6.5.1, 

species cards for 0DK 

• 

REACTIONS 

6.5.2, 

reaction cards for 0DK 


LAST REAX 

INESTS 6. 3. 2. 6, inerts cards for |Z^. (to inert those species that 

. are not named in reactions) 
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Card Input 

Section No., 

Description 

THIRD BODY 

• 

• 

LAST CARD 

6.5.2.7, 

reaction rate ratios for third holies 

$ODK 

SEND 

6 . 5.3 

0[K module namelist 

STRANS 

SEND 

6.6 

TRANS module namelist 

SMOC 

SEND 

6.7 

^t2)C module namelist 

SBLM 

SEND 

6.8 

SJf module namelist 


•feble 6-1 can be used as a guide when preparing Input for given problem. It 
lists the data sets in the order in which they must appear in the data deck, and 
also shows the special cards which must appear in each set (first card, last card, 
etc.) if the program is to function prqperly. The table is basically self-explana- 
tory when used together with the detailed ir^ut descriptions which follow. 

Certain special options to the computer program are described separately In 
Section 6 . 9 . 

An ir^)ut data card listing for a sample case i3 presented in Section 7, follo- 
wed by the corresponding computer output. In preparing input to the computer prog- 
ram it is useful to review this input card listing. 

Successive cases can be run using the computer program but complete data should 
be input for each case. 
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6.1 Ihenaodynamic Data . 

> ( 

Ordlmrlly, a thenoodynaiaic data file Is available for use with the program, 

! and is assigned to loglcial unit 25. Ihe input described here can be used to gene- 

I rate a thermodj'namic data file if one is not available. 

This da^ set is identical to the THE3M) DATA described in i^pendix D of NASA 
SP-273 (i.e. Reference 9). 

Using this data set, thermodynamic data curve fit coefficients may be read from 
cards. The curve fit coefficients are generated by the PAC computer program desc- 
ribed in NASA TN D-4097 (i.e. Reference 22). 


The thermodynamic data (i.e. , etc.) are expressed as functions of tempera- 

ture using 5 least squares curve fit coefficient (a^.j) and two integration 
constants (ag^y) as follows; 


Pip 

IT 




agT a,T^ a a^T^ a, 
“ra ’ * T" * 3 T“ ^ ~ * T" 


S°p a,T^ a.T^ a^T^ 

TT ’ a^ InT + a 2 T + + -y— + -j— + 


I 

For each species, two sets of coeffici«it3 (a-j_y and a^_y) are specified for 
two adjacent temperature intervals, lower and upper respectively. For the data 
available in Reference 9 the lower temperature interval is yxP to 1000°K and the 
upper tanperature interval is 1COO°K to 5000°K. 


Ref. 22. McBride, B.J., and Ciordon, S. , "Fortran IV Program for Calculation of 
Thermodynamic Data", NASA TN-D-4097, Aug. 1967. 
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The input format required for this thermodynamic data is defined in Thble 6-2. 
Data cards for the species AR, H, H 2 » H 2 O, N 2 » 0, OH, and O 2 are listed in Thble 
6-3 as examples. Thermodynamic data coefficients for mar^ chemical species are 
supplied with the computer pr gram, a list of these species is presented in Table 

6-4. 


Data Thpe Genera ticxi and Usage: 

A computer run using thermodynamic data card irput will generate a data ta^ on 
logical unit JANAF. This tape may then be saved and used at a later time. The 
program writes the THE2M) data card Images on unit JANAF as read but with 1»o minor 
exceptions. The THERMO code card and the card numbers in card column 80 are 
omitted. 


If tdiermodynamlc data cards are not iiput, the program assumes the thermo- 
dynamic data is on logical unit JANAF. Logical unit JANAF is curr^tly assigned a 
value of 25 • 
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TABLE 6-2 FORMAT FOR THERMODYNAMIC DATA CARDS 


Card 

order 

Contents 

Format 

Card 

column 

1 

THERMO 

3A4 

1 to 6 

2 

Temperature ranges for 2 sets of co- . . 
efficients: lowest T, common T, and 
highest T 

3F10.3 

1 to 30 

3 

Species name 

3A4 

1 to 12 


Date 

2A3 

19 to 24 


Atomic symbols and formula 

4(A2,F3.0) 

25 to 44 


Phase of species (S,L, or G for solid, 
liquid, or gas, respectively 

A1 

45 


Temperature range 

2F10.3 

46 to 65 


Integer 1 

ns 

80 

4 

Coefficients a.* (1 ° I to 5} in equations 
(6-1) to (6-3) ‘ 

5(E15.8) 

1 to 75 


(for upper temperature Interval) 
Integer 2 

15 

80 

5 

Coefficients in equations (6-1) to 6-3) 
(al, al for upper temperature interval 
ana a^, a 2 , and a^ for lower) 

3(E13.8> 

1 LU /O 


Integer 3 

IS 

80 

6 

Coefficients in equations (6-1) to (6-3) 
(a.,ae,ae,a_ for lower temperature in- 

te^an ® ^ 

4(E15.8) 

1 to 60 


Integer 4 

120 

80 

(a) 

Repeat cards numbered I to 4 in cc 80 
for each species 



{Final 1 
card) 

END (Indicates end of thermodynamic 
data) 

3A4 

1 to 3 


^Gaseous species and condensed species with ohiy one condensed phase can be 
in any order. However, the sets for two or more condensed phases of the same 
species must be adjacent. If there are more than two condensed phases of a 
species thei: sets must be either in increasing or decreasing order according 
to their temperature intervals. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 

TABLE 6-3. THERMO DATA CARDS FOR AN O2/H2 PROPELLANT 
(Species AR, H, H2. H2O, N2, O, OH. and O^) 


THERMO 

300«000 1000*000 5000*000 

AR L 5/66AR 100 000 000 06 300*000 5000*000 


0.35000000C 01 0* Q* 0* 0* 

•0*74537502E 03 0,436b0006E 01 o*250D0O0OE 01 0* 0* 

0* 0* •0*7453749RE 03 0.43660006E 0] 

H J 9/65H 100 000 000 06 300*000 5000*000 

o* 2 Sooooooc 01 0* 0* 0* 0* 

0*?S471627E 05»0*4601 1763E 00 n*2S000000E 01 0* 0* 

0* 0* 0.25471627E 05-0*4b0ll762E 00 

H2 J 3 / 6 IH 20 00 00 06 300*000 5000*000 

0*31001901E 01 0*511l9464E-03 0*52644210E-07-0*34909973E-10 0*36945345E-14 
•0*87738o42E 03»0* 19629421E 01 o*3057445iE 01 0*2b765200E»02-0*58099i62E-05 
O*552lO39lE-O0-O*18122739E-ll-ft.98B9O474E 03-0*22997056E 01 
H20 J 3^61H 20 100 000 06 300.000 5000.000 

C.27167633E 01 0*29451374E*02-0*80224374E«06 0*10226682E*09«0,48472145E*14 
-0.2990SR26E 05 0.66305671E 01 n*40701275E 01-0. 1 1084499E«0? 0*41521 180E»05 
:\29637404E-08 0,90702103E-12-ft,30279722E 05*0,32270046E 00 
N2 J 9/6SN 20 00 00 OG 300*000 5000*000 

0*28963194E 01 0, 15154866E-02-0.57235277E-06 O*990O7393E-1O-O*6S2235S5E-14 
•0*90586184E 03 0.61615148E 01 0.36748261E 01-0* 12091500E*02 0* 23240102E-U5 
•0*63217559E-09-0*22577253E-12-o*10611588E 04 0.23580424E 01 
0 J 6/620 100 000 000 06 300*000 5000*000 

0*2542o 596£ 01-0.27550619E..04-0.31028033E-08 0.4551 0674E«U.0*436805l5E-15 
0.29230803E 05 0*49203080E 01 n*29464287E 01-0* 16381665E-02 0*24210316E*05 
•0*16028432E*08 0*38906964E-12 0*29]47644E 05 0*29639949E 01 
OH J 3/660 IH 100 000 06 300*000 5000*000 

0*29106427E 01 0,95931650E-n3-o,l944l702E.06 0. 1 3756646E-10 0* l4224542E*15 
0*39353815E 04 0,54423445E 01 0*38375943E 01->0* l0778856E-02 0*96630378E«06 
0*187l3972E-09-0*2257l094E-12 0.36412823E 04 0.49370009E 00 
02 J 9/650 20 00 00 06 300*000 5000*000 

0.36219535E 01 0,736l8264E-03-0*l9652228E-06 0*36201558E-10-0.28945627E*14 
•0*12019825E 04 0*36150960 e 01 0*362559BSE 01.0.18782l84E<02 0* 70554544E-05 
-O,67635137E-O0 0,21555993 e«1 1«o, 10475226E 04 0.43052778E 01 
fNO 
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Table 6-4. SPECIES WITH THERMODYNAMIC DATA PROVIDED 


AL(S) 

BCL 

BEO(S) 

AL(L) 

BCL+ 

BEO(L) 

AL 

BGLF 

BEO 

AL+ 

BCL2 

BEOH 

ALB02 

BCL2'(’ 

BEOH+ 

ALCL 

BCL2- 

BE02H2 

ALCL+ 

BCL3 

BE20 

ALCLF 

BF 

BE20F2 

ALCLF2 

BF2 

BE202 

ALCL2 

BF2+ 

BE303 

ALCL2+ 

•BF2- 

BE404 

ALCL2- 

BF3 

BR 

ALCL2F 

BH 

BR2a) 

ALCL3(S) 

BHF2 

BR2 

ALCL3(L) 

BH2 

C(S) 

ALCL3 

BH3 

C 

ALF 

BN(S) 

C+ 

ALF+ 

BN 

C- 

ALF2 

BO 

CCL 

ALF2+ 

BOCL 

CCL2 

ALF2- 

BOF 

CCL3 

ALF3(S) 

B0F2 

CCL4 

ALF3(S) 

B02 

CF 

ALF3 

B02- 

CF2 

ALH 

BS 

CF3 

ALN(S) 

B2 

CF4 

ALN 

B20 

CH 

ALO 

B202 

CH2 

ALO+ 

B203(L) 

CH20 

ALOCL 

B203 

CH3 

ALOF 

B303CL3 

CH4 

ALOH 

B303F3 

CN 

ALOH+ 

BE(S) 

CN+ 

ALOH- 

BE(L) 

CN- 

AL02 

BE 

CN2 

AL02- 

BE-l- 

CO 

AL02H 

BEB02 

CCCL 

AL2CL6 

BECL 

CCCL2 

AL2F6 

BECL+ 

CCF 

AL20 

BECLF 

COF2 

AL20+ 

BECL2(S) 

COS 

AL202 

BECL2(L) 

C02 

AL202+ 

BECL2 

C02- 

AL203(S} 

BEF 

CP 

AL203(L) 

BEF2(S) 

CS 

AR 

BEF2(S) 

CS2- 

AR+ 

BEF2a) 

C2 

B(^ 

BEF2 

C2- 

ea) 

BEH 

C2CL2 

B 

BEH+ 

C2F2 

B+ 

BEN 

C2F4 


C2H 

FECL2(S) 

KOH 

C2HF 

FECL2IL) 

vrtufe\ 

C2H2 

FECL2 

KOHCL) 

C2H4 

FECL3(S) 

K2 

C2N 

FECL3(L) 

K20(S) 

C2N2 

FECL3 

LI(S) 

C20 

FEO(S, 

LI(L) 

C3 

FEOCL) 

U 

C302 

FEO 

LI+ 

C4 

FE02H2(S) 

UCL(S) 

C5 

r'E02H2 

UCLa) 

CL 

FE03H3(S) 

LICL 

CL^ 

FE203(S) 

LIF(S) 

CL- 

FE304(S) 

LIF(L) 

CLCN 

H 

LIF 

CLF 

H+ 

LIF2- 

CLF3 

H- 

LIFO 

CLO 

HALO 

LIH(S) 

CL02 

HBO 

LIH(L) 

CL2 

HBO+ 

LIH 

CL20 

HB02 

LIN 

CS(S) 

HCL 

LIO 

CS(L) 

HCN 

LIO- 

CS 

HCO 

UOH(S) 

CS+ 

HCO+ 

LIOHCL) 

CRCT.ISl 

HOP 

LIOH 

CSCLCB) 

HF 

LION 

CSCLa) 

HNO 

U2 

CSCL 

H02 

U2CL2 

CSF(S) 

H2 

LI2F2 

CSF(L) 

H20(S) 

LI20{S) 

CSF 

H20a) 

LI20(L) 

CSO 

H20 

U20 

CS2 

H202 

U202 

CS2CL2 

H2S 

LI202H2 

CS2F2 

H3B306 

U3 CL3 

CS2C 

HE 

LI3F3 

E 

HE+ 

MG(S) 

F 

K(S) 

MG(L) 

F- 

KCL) 

MG 

FCN 

K 

MG+ 

FO 

K+ 

MGCL 

F02 

KCL(S) 

MGCL+ 

F2 

KCL CL) 

MGCLF 

F20 

KCL 

MGCL2(S) 

FE(S) 

KF(S) 

MGCL2(L) 

FE(S) 

KF(L) 

MGCL2 

FE(S) 

KF 

MGF 

FE(L) 

KF2- 

MGF2(S) 

FE 

K2F2 

MGF2(L) 

FECL 

KO 

MGF2 



Table 6-4. (cont'd) 



MGH 

0 

SIN 

MGN 

0+ 

SIO 

MQO(S) 

0- 

SI02(S) 

MGO(L) 

OH 

SI02(S) 

MGO 

0H+ 

SI02(S) 

MGOH 

OH- 

SI02CL) 

MGC/H+ 

02 

SI02 

MG02H2 

. 02- 

SIS 

N 

P 

SI2 

NF 

P(S) 

SI2C 

NF2 

P+ 

SI2N 

NF3 

PCL3 

SI3 

NH 

PF3 

XE 

NH2 

PFS 


NH3 

PH 


NO 

PH3 


NO+ 

PN 


NOCL 

PO 


NOF 

PS 


NOF3 

P2 


N02 

P4 


N02- 

S(S) 


N02CL 

SCL) 


N02F 

S 


N2 

s-i- 


N2C 

SF4 


N2H4 

SF6 


N20 

SH 


N204 

SN 


NA{S) 

SO 


NACL) 

SOF2 


NA 

S02 


NA+ 

S02F2 


NACL(S) 

S03 


NACL(L) 

S2 


NACL 

SI(S) 


NAF(S) 

sia) 


NAF(L) 

SI 


NAF 

SI+ 


NAF2- 

SIC 


NAH 

SIC2 


NAO 

SICL 


NAO- 

SI0L2 


NAOH(S) 

SICL3 


NAOH(L) 

SICL4 


NAOH 

SIF 


NA2 

SIF2 


NA2CL2 

SIF3 


NA2F2 

SIF4 


NE 

SIH 


NE+ 

SIH4 
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6.1.1 THERMODYNAMIC DATA BELOW 300°K. 


Ordinarily this data set is not required. However, for low temperature calcu- 
iatlons it iBay be necessary to extend the curve fit data in the Thermodynamic Data 
file (see Sectic« 6.1). The lower temperature limit, T^, in the Thermodynamic Data 
supplied with the program is 300^K. 

Thermodynamic Data below the temperature, T^, may be input by data cards as 
described below. 


card 1 L0W T CPHS 

card 2 n 


0 0 0 0 

card 3 T^K Cp^ Hj, Sj 1 


Directive for start of low ten^)erature 
CPHS tables (col 1 throu^ 10). 

12 character species name, left Justi- 
fied, followed 1:^ the integer, n, 
pinched in colunn 21. The itegrer n 
must be such that 1<n<3 and represent 
the nunter of Thermod^nmic Data 
poLiuc to be input for this species. 


First Thermodynamic Data point for the 
above species, irput 4F 10.0, 15. 


card n+2 T?K Cp hS s2 nth (Kn<3) nth Thermodynamic Data point for the 

T ^ above species, input 4F10.0, 15. 

. . . Repeat cards through n+2 above for each species to be irput. 

Tanperature must be T <T <T <T . 

^ Z 3 I 

(final ERD L0W T CPHS end directive (col 1 through 14) 

card) 
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An example of this irput is given in Table 6-5 which shows a card listing extending 
the Thermodynamic Data for an propellant to 100°K. Data in Tible 6-5 is 
taken directly from JANAF tables (Reference 23) > except for Argon which is 
taken from NASA SP-3001 . 


The quantity is defined as 


hJ ■ - hJ , cal/mole 


298 


and 


c2 , cal/mole - deg K 


, cal/mole - deg K 


Ref. 23. Stull, D.R. , Prophet, H. , et al., JANAF Theimochemical Thbles, Second 
Edition, NSRDS-NK National Standard Reference Data Series, National 
Bureau of Standards, June 1971. 


* j 
" i 
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TABLE 6-5. LOW lEMPERATURE DATA FOR AN O 2 /H 2 PROPELLANT 


> 

LOW T CPHs 





&R 


2 




100,0 

4,9681 

•984.5 

31*556 

1 


200,0 

4,9681 

-487.7 

34,999 

2 


H 


2 




100.0 

4.968 

51118. 

21,965 

1 


200,0 

4,968 

51614. 

25.408 

2 


H 2 


2 




100.0 

5.393 

•1265. 

24.387 

1 

> 

200.0 

6,518 

- H 2.0 

28.520 

2 


h20 


2 




100.0 

7.961 

•59378,9 

36.396 

1 

* 

200.0 

7,969 

•58581,9 

4 i ,9 l 6 

2 


N 2 


2 




100.0 

7,074 

•1387.0 

38.113 

1 


200,0 

6.989 

-684, 

42,986 

2 


0 


2 



*• 

100.0 

5.666 

58479. 

32.466 

1 


200.0 

5,434 

59036, 

36.340 

2 


OH 


2 




100.0 

7,567 

7879, 

35.852 

1 


200,0 

7,309 

8623. 

41*021 

2 


02 

' 

2 




100.0 

6.958 

•1381. 

41*395 

1 


200*0 

6*961 

-685. 

46*218 

2 

’ 

end uo « t 

cPhs 
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6.2 


Title Cards 


This Input permits labeling of runs with alphanumeric information. As mar^ 
title cards as desired may be input in sequence. Card format is as follows: 

col 1-5 col 6-77 

TITLE aiv alphanumeric information 

It is not necessary to input title cards. 
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6.3 DATA Directive and $DAIA Namelist Input . 

Ihe DATA directive and the $DATA Namelist input set. described below must always 
be input. It is required for all problems since it contains the iiput that cont- 
rols which calculation modules are to be executed. 

The first liput item must be a single card, called the DATA directive card. 
The format of this card is as follows; the letters DATA must be pundied in coluims 
1 throu^ 4 . The DATA card is used to inform the program that the $DATA namelist 
input is to follow. 

Ihe card following the DATA card must contain the name $OATA, and all cards in 
the namelist irput set must start in column 2 or greater. Since Namelist input is 
card interpretive, items can be input in any order. last card in the set must 
contain SEND. 

Users unfamiliar with Namelist input are referred to their P0FO31AN reference 
manual. 


Assumed 

Item Description Ifaits Value(s) 

DATA DATA directive card 

$DATA Namelist name, read in Subroutine 

PR0BLM 
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6.3.1 Specification of Modules to be Executed . 


If a module is to be executed, it is necessary to indicate the fact by input of 
a module flag as described below. For example, if a problem requires that the 0DE 
tsodulo be ran, it is necessary to input 0Lfi * 1. Oily certain combinations of 
modules are allowed. These are described in Table 6-6. 'Ihe module flags are: 


Assumed 


Item 


Description 

Ihlts 

Value(s) 

0D& 

3 

Set 0DE a 1 if the 0£® module is to be 
executed. 

none 

0. 

0OK 

3 

Set 0DK - 1 if the ODK module is to be 
executed. 

none 

0. 

TE® 

S 

Set IDE « 1 if the TEE option of the 
module is to be executed 

none 

0. 

TDK 

3 

Set TDK ■ 1 if the TDK option of the 
V$C module is to be executed 

none 

0. 

TDP 

m 

Set TDF ■ 1 if the IDF option of the 
VHtC is to be executed 

none 

0. 

TDKIL 

a 

Set TDKIL * 1 if the TH( option of the 
\/!$C is to be executed with an initial 
data line input as described in 
Section 6.9.2. 

none 

0. 

BLM 

3 

Set BLM . 1 if the BLM option of the 
M0C module is to be executed 

none 

0 . 

PPG50PT 

3 

Set PFD0PT a 1 if the constant proper- 

none 

0. 


ties option of the MOC module is to be 
execute. See Section 6.9.1 for input 
instructions. 
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•Dable 6-6: Usage of the Module Flags 


Input 


Mode of Execution 


0DE-1, 


0DE-1, 0DK*1, 


0DK«1, 


0DE-1, 0DK«1, TDK»1, 


0DE-1, 0DK-1, TDK-1, BLM-1, 

0DE-1, 0DK-1, TDK-1, BLM-1, 
IRPEAI-1, 


0DE-1, 0DK-1, TDK-1, BLM-1, 
IRPEAT-2, 

TDKIb-1, 


TDE-1, 

lDE-1, BLM-1, 

TOE-1, BLM-1, IRPEAT-1 or 2, 


0I£ is run alone. Options other than the 
tTScket (HKi'-T) option are allowed. See 

Section 6.4 >3 and Reference 9, 

0DE is run. 0DK Is run with 0DE providing 
start conditions. See Section 6.5.1. 

0DK is run alaie with initial conditions 
supplied by the user. See Sections 6.5.1 and 

6.5. 1.1. 

is run. 0DK i„ run with 0DE providing 
start conditions. TDK is run wilh 0DK and 
TRANS providing M0C initial data line 

conditions. The number of 0DE and 0DK runs 
will be equal to NZ0NES. 

As above with a BUI run added. 

As above with 0DE, 0DK, and TDK repeated 
using the BLM results. See Section 6.8.4. 
This option is for including the BLM resulte 
with regen cooling. 

As above except that the 0DE and 0DK runs are 
not repeated because there is no heat input 
correction due to regen cooling. 

TDK is run with a MfSC initial data line 
input as described in Section 6.9.2. 

0DE will be run for NZ0NES. TDE is run with 
0TR chemical equilibrium gas properties, and 
with TRANS providing M0C initial data lini 
conditions. 

As above with BU!. 

As above with TDE repeated. 
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Ifeible 6-6: Usag-s of the Module Flags (continued) 




Input 

0DE*1, 0DKa1, TDF»1, 

0DE»1, 0DK-1, TDF=1, 

0DE»1, 0DK»1, TDF«1, BLM-1, 
IRPEAT-1 or 2, 

BLM-1, 

PPG0PT-1, 


Mode of Execution 

0[£ is run. 0DK is run with 0E£ providing 
start coTiditiaio. However, the 0DK run will 
have a frozen chemical composition. TDK wiJ 1 
also be run with a frozen composition. 

As above with HM. 

As above with TDF repeated. 


BLM is run alone. This option requires a 
large amount of input that is provided 
automatically when BLM is run after TIE, TDK 
or TDF. 

The ILIANS and modules are run with 

constant gas prc^rties. See Section 6.9*1. 




6.3.2 Inputs for Control for -ttie Program . 


Item 

IRPEAT 


NZ0NES 


SI 


I0FF 


Description Units 

Set IRreAT » 1 or 2 to request that a none 

IKC (or TDE) calculation be automati- 
cally repeated after the BLU module 
has been used to calculate a displaced 
nozcle wall. 

If IRPEAT » 1 , the (9DE and 0CK module 
executions will be repeated with 
adjusted enthalpies for regen cooling, 
see Section 6.8.4. 

If REPEAT 2, the 0DB and 0DK module 
executions will not be repeated. !Die 
system enthrJpy will be unchanged. 

Nuipber of zones, N, to be used In the none 

anlysls. The and iXK modules will 
be executed N times. Zone 1 repre- 
sents the flow adjacent to Ihe nozzle 
axis, and Zone N r^resents flow adja- 
cent to -ttie nozzle wall. 

If SI * 0, Ehgllsh units are to be none 

used for Input and output. If SI > 1 , 

SI units are to be used for Irput and 
output. The SI units required for 
Input are shown In parenthesis. 

Vftven BLM Is run with liie MX module, 
values deflnli;g the boundary layer 
edge conditions; x@» Yq* Uq, Tg, and 
Pg are autcmatlcally calculated and 

stored In the XIN0, RIN0, UE0, TB0, 
and PB(8 arrays. If ICFP Is not Input, 
then these values are stored starting 
with the first entry In each array. 

If liZ^ Is Input then these values are 
stored starting In the I0FF + 1 entry 
of each array. Values must then be 
Input Into entries 1 through I(3FF. 

The boundary layer calculations will 
start at position XIN0(1) of these 
arrays using the $BIM Irput. This 
irpul allows the user to account for 
the development of the boundary layer 
in the nozzle changer, upstream of 
BC31AT. I0FF < IX. 

If I0FF is irput, do not input NXIN0 
in Bin. 


Assumed 
Value (s) 

0 


1 


0 

0 
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6.3*3 Specification of nozzle geometry . 

To eliminate redundancy, all geometry iiputs that are common throu^ut the 
modules are ir^ut here. Note that the center of the nozzle geometry coordinate 
system ia at the centerlLne at the throat plane, and that all coordinates are 
normalized by the nozzle throat radius, r^. hence, axial positions upstrerm of the 
throat are always negative numbers. In the figures describing the geometry, posi- 
tive angles ure shown as counter-clockwise, and negative angles are shown as clock- 
wise. 


Geometric area ratios at which 0DE and 0DK print out is to be made are 
specified using the iiput arrays ASUB(1) and ASUP(1). 

Ite nozzle geometry is defined in Figure 6-1. Hie calculations start at 
the downstream end of tl^ combustion chamber with a subsonic area ratio of ECRAT, 
as shown. The circular arcs RI and RWTU cannot overlap. Thus, it is necessary 
that ECRAT, RI, RWHJ, and ThETAI be input such That 

v/ ECRAT > 1 + (RI + RWTU) (1 - cos THEJTAI). 


If this condition is not met, subroutine PRES of 0DK will print the terminal error 
message: 


INLET GEOMETRY INCOMPATIBLE WITH INITIAL CONDITIONS. 

In addition, the trarisonic analysis requires that a value of RWTU _> .5 be 
Ir^ut. 

The wall geometry downstream of the nozzle throat can be specified using any 
one of several options. All of these geometries begin with a circular arc of ra- 
dius RWTD which extends from the throat point through an angle of THETA. Geomet- 
ries that can be iiput to both the 0DK and M0C modules are described in Secticxi 

6.3.3.I. 
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(0,0) (0,xJ 


Figure 6-1: Nozzle Geometry, all coordinate 

values are normalized by RSI = r^. 
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0 


!0 


Item 


Description 

Units 

Assumed 

Value(s) 

RSI 


Nozzle throat rad is is, r^ 

in, 

(meters) 

0. 

ASUB(I) 

S 

Subsonic area ratios at ^ch informa- 
tion will be printed. 

none 

0. 

NASUB 

= 

Number of entries in the ASUB array 
< 50. 

none 

0. 

ASUPd) 

SB 

Supersonic area ratios at which infor- 
mation will be printed. For the TCE 
option, ASUP aaist contain an entry 
the nozzle expanslcxi ratio. 

none 

0. 

NASUP 

SB 

I*junber of entries in ttie ASUP array 
< 50. 

none 

0. 

BCRAT 

a: 

Nozzle inlet contraction ratio for use 
in 01E and 0DK calculations. 

none 

0. 

RI 

ss 

Normalized inlet wall radius. 

none* 

0. 

THETAI 

at 

Nozzle inlet angle. 

degrees 

0. 

RWTU 

» 

Upstream normalized wall throat radius 
IWTO__> .5 is required. 

none* 

0. 

rWPE 

SK 

Hype of nozzle wall to be itput. 

none 

0. 


ITyPE a 0, if the real wall contour is 
irput. 


ITfPE » 1, if the potential flow wall 
contour is iiput. 

If niPEAT a 1 or 2, then the nozzle 
wall will be displaced by +fi* as cal- 
culated by BLM when the TESf (or TIE) 
calculations are repeated. This dis- 
placement is -5* for the real wall 
(ITYPE * 0), or +6* for the potential 
flow wall (ITyPE a 1). 


* Normalized by the throat radius, r^ 


! 

A 
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I 


6 . 3 . 3.1 

Exhaust Nozzle Geometry Specification, for 0OK and the M|2)C Modules. 

» 

Item 


Description 

Units 

Assumed 
Value (s 

IWALL 

ret 

option flag for specifying the 
downstream wall. 

none 

0 



1 cone (^tion (input RVID, THETA, and 
EPS) 




= 

2 parabolic nozzle contour c»tion 
(input RWTD, THETA, RMAX, 2MAX) 




s 

3 circular arc nozzle contour option 
(input as for IWALLs2) 




a 

4 nozzle contour (spline) qjtion (in- 
put ROTD, THETA, THE, RS, ZS, NWS) 




a 

5 cone with specified end point, 
(input RWTD, BMAX, and ZMAX) 




a 

6 ideewed parabola option. (ir^ut 
RWD, THEIA, RMAX, ZMAX, THE) 



The items required for the various IWALL options are: 



Rwro 

a 

downstream wall throat radius of cur- 
vature ratio** 

none* 

0. 

THEIA 

3 

nozzle attachment angle 



THE 

3 

nozzle exit angle (input if IWALLs4, 
or 6) 

deg. 

0. 

EPS 

a 

nozzle expansion ratio (iiput if 
IWALU1 only) 

deg. 

0. 

RMAX 

» 

normalized radius at the nozzle exit 
plane (input if IWALL»2 or 3) 

none* 

0. 

ZMAX 

a 

normalized axial position at the 
noiLzle exit plane (input if IWALL » 
2 or 3) 

none* 

0. 


* Normalized by the throat radius r^ 

** If a comer expansion (i.e, Prandtl-Meyer fan) is desired, a value of RWTD « 
.05 is recommended. Experience has shown that values smaller than this give the 
same result but are computationally less efficient. 


6-22 



V. 


i 

i 



K 


Assumed 


Item 

Description 

tJhits 

Value(s) 

RS(2) 

table of normalized wall radii down- 
stream of the nozzle tangency point 
(irput if riiALL«4)« Tiie input tables 
RS and ZS start with the second entry 
because the first entry is calculated 
autcuatically by the program, l.e., 
RS(1)»r,j, and T ZS(1)»z,j, in Figure 6-2. 

none**’ 

0. 


Ihe T wall angle at this position is 
also calculated so that the spline 
contour will be properly joined to the 
nozzle throat contour. 



ZS(2) 

table of normalized axial position 
downstream of the nozzle tangency 
point (input if IWALL«4)< 

none**’ 

0. 

NWS 

total number of entries in the RS, ZS 
tables. Includes the first entry NWS 
< 50 (input if IWALL=4). 

none 

0 

RZN0RM 

Optional normalizing factor for the 

none 

1. 


ZSy RlNy ZIN tables. For exan^le, 
if RSy ZS, and RIN, ZIN were input as 
dimensional numbers, RZN0RM would be 
the throat radius in those units. 


See the input variable RZNI2KV1. 
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6.4 GfDE INPUT DATA (ALL PROBLEMS SPECIFYING QfPE) 


The CfDE Input data described here Is exactly as defined In NASA SP-273, 
Reference 9, except namelists input $INPT2 and $RX1NP have been combined into 
a single list named $0DE. Any type of equilibrium calculation available with the 
computer program described in Reference 9 can thus be computed using the $(2TDE 
Input data*. In this document, however, only the RKT option of namelist is 
described. The RKT option differs from that of Reference 9 for problem types 
other then single zone 0DE. 

The QfDE input data consists of the following input groups: 

1. REACTANTS directive card, followed by up to 15 

data cards, followed by a blank card, 
specifying reactants. 

2. QfhAlT and INSERT directives to omit or insert species for 

equilibrium/frozen calculations . 

3. NAMELISTS directive card followed by input 

namelist $(2^E specifying input 
case data . 

6.4.1 REACTANTS CARDS 


This set of cards is required for ail problems. The first card in the 
set contains the word REACTANTS punched in card columns 1 to 9. The last card 
in the set is blank. In between the first and last cards may be any number of cards 
up to a maximum of 15, one for each reactant species being considered. The 
cards for each reactant must give the chemical formula and the relative amount 
of the reactant. For some problems , enthalpy values are required. The format 
and contents of the cards are summarized in Table 6-7. A list of some REACTANTS 
cards is given in Table 6- 8 

Relative amounts of reactants. - The relative amounts of reactants may be 
specified in several ways. They may be specified in terms of moles, mole fraction, 
or mole percent (by keypunching M in card column 53) or in terms of weight, weight 
fraction, or weight percent (blank in column 53). 

Relative amounts of total fuel to total oxidants can also be input. For this 
situation, each reactant must be specified as a fuel or an oxidizer by keypunch- 
ing an F or O, respectively, in column 72 of the REACTANTS card. The amounts 

♦These options include TP,HP,SP, TV, UV, or SV problems, Chap.nan-Jouguet de- 
tonation problems, and incident or reflected shock problems. 
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Of POOR QUALITY 


TABLE 6-7 REACTANTS CARDS 


Order 

Contents 

Format 

Card columns 

First 

REACTANTS 

3A4 

1 to 9 

Any 

One card for each reactant species 
(maximum 15). Each card contains: 




(1) Atomic srmbots and formula num- 
bers (maximum S sets)^ 

«A2. F7.5) 

1 to 45 


(2) Relative wei^ht^ or number of 
moles 

F7.5 

46 to 52 


(3) Blank if (2) is relative weight or 
M if (2) is number of moles 

Al 

53 


(4) Enthalpy or internal 
cal mole 

F9.5 

54 to 62 


(5) State: S. L* or C for solid, 
liquid or t:as« respectively 

Al 

63 


(6) Temperature associated with 
enthalpy in (4) 

F8.5 

64 to 71 


(7) F if fuel or 0 If oxidant 

Al 

72 


(8) Density in f* cn^^ (optional) 

F8.5 

73 to 63 

Last 

Blank 




^ProKrara will calculate the enlhalpv or intcrnil energy (4) for species in 
the THERMO data at the temperature (6) if zeros are pimched in card 
columns 37 and 38. (See section Reactant enthalpy for additional in* 
furnuition.) 

^eLntive weight of fuel in total fuels or cxidant in total oxidants. All 
rcacUinis must be j^iren either all in relative weights or all in number 
«if moles. 
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OFilcaHiAL E ;Vj'£ Eii 

OF POOR QUALITY 


lASLE 6*8 LISTING OF SAMPLE REACTANTS CAROS 


R6ACTANT8 



H 2. 

too. 

0. 0298,19 f 

N ,7806810 .209795AR, 084682 

too. 

-7.2021840298,19 0 


WUCT8NTS 
Ml. M 

4. CLl, 0 4, 

72.08 

-70730. 

8296,19 

p 

Cl, H 

1.889990 .031298S .008419 

18.58 

•2998, 0821298, 15 

p 

4U. 


8,00 

♦CoO 

S29a,lS 

f 

M61, 0 

1. 

.20 

•143700. 8298,19 

f 

M2. 0 

1. 

• 16 

•68317,4 1298.19 

f 


RK8CT8NTS 
M 7. 

00 . 

lOOo 

0. 

6298,19 

f 


0 2. 

00 

lOOoO 

0.0 

6298.19 

0 


REACT ANT8 

M2, H 6. C 2, 


50.0 

12744,8 

1298.19 

P 

,788 

M2. K 4. 


50.0 

12050. 

1296,19 

P 

1.003 

R 2, 


100. 

•3030.8921 89,24 

0 

1.94 

REACT AMT8 
111. 


100. 

0. 

8298.15 

p 


R 2. 


100. 

•3030.892c 89,24 

0 

1,94 

reactants 

M2, H 4. 


00. 

12100. 

L298,15 

p 

1.003 

BEl, 


20. 

0.0 

8258.13 

p 

1.89 

M2. 0 2. 


100. 

>44880, 

L298.15 

0 

1.407 


^Listad above an six txamplts. Each axampla must and with a blank card, 


given on the REACTANTS cards are relative to total fuel or total oxidant rather than 
total reactant. 

There are four options in the $0DE namelist for Indicating relative 
amounts of total fuel to total oxidant as fcllow>5: 

1 . Oxidant to fuel weight ratio (0F is true) 

2. Equivalence ratio (ERAH0 is true) 

3 . Fuel percent by weight (FPCT is true) 

4. Fuel to air or fuel to oxidant weight ratio ^A is true) 

For each option, except 0DE with NZ0NES«1, the values are given In the 
0FSKED array of $0DE (described in Section 6.4.3). For 0DE with NZ0NES»1, 
the MIX array is used, as described in Reference 9 . 


Reactant enthalpy. Assigned values for the total reactant are calculated 
automatically by the program from the enthalpies of the Individual reactants. 
Values for the individual reactants are either keypunched on the REACTANTS 
cards or calculated from the THERM0 data as follows; 

Enthalpies are taken from the REACTANTS cards unless zeros are punched 
in card columns 37 and 38. For each REACTANTS card with the “00" code, an 
enthalpy will be calculated for the species from the THERM0 data for the tempera- 
ture given in card columns 64 to 71. 

When the program is calculating the individual reactant enthalpy for 
values from the THERM0 data, the following two conditions are required; 

1. The reactant must also be one of the species in the set of THERM0 data. 

For example, NH^(g) is in the set of THERM0 data but NH^(t) is not. Therefore, 
if NH^(g) is used as a reactant its enthalpy could be calculated automatically, 
but that of NH^(i) could not be. 

2. The temperature T must be in the range Tj^^l . 2 s T s T^^^^^x 1.2 where Tj^^ 
to Is the temperature range of the THERM0 data. 

For cases with NZ0NES > 1 (see Problem card, Section 6.3) it may be de- 
sirable to modify the enthalpy of each zone. This can be done by using the DELH 
input array. For the i^ zone the i'^ DELH entry will be added to the system 


enthalpy as computed by 0DE from the reactants cards (see above). For example, 
overall system enthalpy of the propellants In the tank can be input through the reactants 
cards and the work added or extracted per zone can be Input by the DELH entries. 

An alternate method would be to input zero enthalpy on the Reactants cards and 
input enthalpy per zone by the DELH entries. 

6.4.2 0MIT and INSERT Cards 

0MIT and INSERT cards are optional. They contain the luimes of particular 
species in the library of Thermodynamic Data for the specific purposes discussed . 
below. Each card contains the wcxd 0MIT (in card columns 1-4) or INSERT (in 
card columns 1-6) and the names of from 1 to 4 species starting in columns 16, 

31, 46, and 61. The names must be exactly the same as they appear in the 
THERM0 data. 

6.4.2. 1 ^MIT Cards 

These cards list species to be omitted from the THERM0 Ham . If 0MIT 
cards are not used, the program will consider as possible species all those 
species in the THERM0 data which are consistent with the chemical system 
being considered. Occasionally it may be desired to specifically omit one or 
more species from considerations as possible species. This may be accomplished 
by means of 0MIT cards. 

6. 4. 5.2 INSERT Cards 

These cards contain the names of condensed species only. They have 
been included as options for two reasons. 

The first and more important reason for including the INSERT card option is 
that, in rare instances, it is impossible to obtain convergence for assigned enthalpy 
problems (HP or RKT) without the use of an INSERT card. This occurs when, by 
considering gases only, the temperature becomes extremely low. In these cases, 
the use of an INSERT card containing the name of the required condensed species 
can eliminate this kind of convergence difficulty. When this difficulty occurs, 
the following message is printed by the program: "LOW TEMPERATURE IMPLIES 
CONDENSED SPECIES SHOULD HAVE BEEN INCLUDED ON AN INSERT CARD". 


The second and less important reason is that if one knows that one or 
several particular condensed species will be present among the final equilibrium 
compositions for the first assigned point, then a amall amount of computer 
time can be saved by using an INSERT card. Those condensed species whose 
chemical formulas are included on an INSERT card will be considered by the 
program during the initial iterations for the first assigned point. If the INSERT 
card were not used, only gaseous species would be considered during the initial 
iterations. However, after convergence, the program would automatically insert 
the appropriate condensed species and reconverge. Therefore, it usually is 
immaterial whether or not INSERT cards are used. For all other assigned points 
the inclusion of condensed species is handled automatically by the prorrram. 
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6.4.3 SgPE NAMELIST INPUT 


The QfDt subprogram contains namelist Input sections $0D£ 
and $SHKINP. The Namelist $0DE must be preceded by a card with NAMELISTS 
punched In card columns 1-9. 

The $0DE Namelist is required if 0DE =1. or TDK =1, in SDATA as 
described in section 6.3. 

For the 0DE problem type any of nine (9) different equilibrium problems can be 
solved. These are TP,HP,SP,TV,UV,SV,RKT,DETN, and SH(?CK. For the 01DE- 
0DK, 0DE-(2fDK-TDK, or TDE problem type, only the RKT problem can be 
solved. In this section only the RKT input option is discussed. Reference 9 is 
.to be used to prepare input for the other equilibrium problems. 

The variables input by the $<?DE namelist are listed in Table 6-9. 

Additional information about some of these variables follows: 

Measure units. - The nronram assumes the pressure in the P schedule to be in 
units of atmospheres unless either PSIA « true, or SI s true. 

Relative amounts of fuel(s) and oxidizer (s) . - These quantities may be specified by 
assigning 1 to 15 values for either o/f, %F, f/a, or r. If no value is assigned 
for any of these, the program assumes the relative amounts of fuel(s) and oxldizer(s) 
to be those specified on the REACTANTS cards. (See discussion in REACTANTS 
Cards, Section 6.4.1) 

RKT problem . - Only one value for chamber pressure, P, Is to be input for cases 
with N20NES> 1 (see $DATA input, Section 6,3). The stagnation pressure used for 
the 1*^ zone will be the value input for P multiplied by the value input in the 
schedule XP . If not input, all XP entries are assumed equal to one. For TDK type 
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problems, zone one Is taken about the nozzle axis of symmetry and the last 
zone is bounded by the nozzle wall. Similarly, the zone will have a mixture 
ratio equal to the 1^ entry in the Ot^SKED schedule = 

Print out will be given for the chamber pressure condition (l.e. stagnation) 
and the throat condition. Print out may be requested at other conditions by use 
of the PCP schedule and the SUBAR and SUPAR schedules. 

The program will calculate both equilibrium and frozen performance unless 
FR0Z * F or EQL » F are input. If FR0Z = F, only equilibrium performance will 
be calculated. If EQL * T, only frozen performance will be calculated. 



ORIGINAL FAGt fS 
OF POOR QUALmr 


TABLE 6-9 VARIABLES IN $0DE NAMELIST 


Variable 

Nc. of 
entries 

Type 

Value 

befc.e 

read 

Definition and comments 

* 

RKT 

1. 

L 

False 

Rocket problem^ 

P 

26 

R 

0 

Assigned pressures: stagnation pres- 
sures for rocket problems: values in 
atm ua’'.SG PSIA, or SI « ,T. , (see be- 
low) 

SI 

1 

I* 

False 

^Values in P array are in N/m* 

PSIA 

1 

L 

False 

^Values In P array are in psia units 

XP 

50 

R 

1. 

Multipliers for the zone stagna- 

tion pressure (zone 1 - inner zone) 

0F . 


L 

False 

Oxidant to fuel weight ratios are to 
be input® 

ERAH0 

1 

L 

False 

Equivalence ratios are to be l^ put® 

FPCT 

1 

L 

False 

Percent fuel by weight are to be input® 

FA 

1 

L 

False 

Fuel to air weight ratios are to be input® 

0ISKED 

50 

R 

0 

For a Rocket problem , and NZ0NES 1 . 
0FSKED will be used rather than MIX 
(see Reference 9). Relative amounts of 
total oxidant to total fuel are input as 
defined by 0F, ERATI0, FPCT, or FA. 
For 0DE-0DK-TDK and TDE problem 
types these valuec define the oxidant to 
fuel ratios for each zone (zone 1 = 
inner zone) 

DELH 

i 

50 

R 

0 

Corresponding to each zone this value 
will be added to the f ystem enthalpy 
input thru the reactants cards . Units 
are BTU/# if PSIA=,T, , Joule/)ciiogram 
if SI*,T. , otherwise cal/gram. 

DELHI 

50 

! 

R 

0 

i 

1 

1 

Corresponding to each zone this value 
will be added to the system enthalpy. 
These values can be used as a 1st est- 
imate for the heat returned to the main 
combustion chamber by regen cooling 
circuits (CfFC input in $BLM). The BLM 
will recalculate these values and, if 
IRPEAT =i in $DATA, rerun the 0OE, 0OK, 
TDK (or TDE) analysis. Same units as 
DELH. 


*If variable Is set to be trie. 
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Table 6-9 (cont’d) 


1 Variable 

No, of 
entries 

Type 

Value 

bef''re 

;ead 

Definition and comments 

I0NS 

1 

L 1 

False 

Consider ionic species^ 

WFL0W 

1 

R 

0 

Input nozzle mass flow option for C?DE- 
^iDK-TDK or TDE problems. If ? value 
KV 1 -L 0 W is input an expansion with 
this mass flow will be computed. Ihe 
values input for P and XP are used as 
estimates for computing stagnation pres- 
sure for each zone, Ihe program will 
adjust these stagnation pressures to 
obtain the desired nozzle mass flow 
within a tolerance of RELERR, Units 
are Ibs/sec if PSI^,!, otherwise kilo- 
grams/sec. 

RELERR 

1 

R 

,0005 

Relative difference between requested 
and computed mass flow rate, Ihe pro- 
gram stops if this error is exceeded. 

PCP 

50 

R 

0 

Compute and print solutions at these 
values of chamber pressure to pressure 
ratio (entries must be >1,) 

SXJBAR 

50 

R 

0 

Compute and print solutions at these 
values of subsonic area ratios (en- 
tries must 7 ^ 1,) 

SUPAR 

50 

R 

0 

Compute and print solutions at these 
values of supersonic area ratio (en- 
tries must 1) 

ECRAT 

1 

R 

0 

Subsonic area ratio to start QfDK. cal- 
culations with computed equilibrium 
conditions. The SUBAR input table 
must Include an entry equal to ECRAT. 

EQTHST 

1 

L 

False 

To start 0DK calculations with com- 
puted equilibrium conditions at the 
nozzle throat,® 

EQL 

1 

L 

Tnie 

Calculate rocket performance assum- 
ing equilibrium composition during 
expansion^. 

PR0Z 

1 

L 

Thje 

Calculate rocket performance assum- 
ing frozen composition during expan- 
slonb. 

LISTSP 

1 

L 

False 

List names and dates of all species ^ 
residing on thermodynamic data used . 

KASE 

1 

I 

0 

Optional assigned number associated 
with case. 


*if variable Is set to be true, 

^Set variable false if these calculations are not desired, 
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6.4.3. 1 OPTION TO PUNCH TABLES FOR BOUNDARY LAYER PROGRAM INPUT 

(DOES NOT APPLY TO BLM) 

Conditions computed along the nozzle wall can be output as puriched cards 
for Input to the BLIMP, TBL, or MABL boundary layer analysis computer programs. 
These conditions are taken by the boundary layer computer program as being the 
Invlscld flow condition at thw edge of the boundary layer. Tables to be punched 
are: x, y, and P/P_ (i.e. the nozzle wall coordinates and the ratio of pressure 
to chamber pressure along the wall), The tables are punched in NAMELIST for- 
mat readable by BUMP (see the punched card listing given at the end of the sample 
output. Section 7). 

A maximum of 50 entries upstream of the throat are saved and punched. 

The wall point at the end of every characteristic is punched up to a maximum of 
500 total table entries . The user may specify a number by which the punched table 
will be offset. Thus, the first point may be output with identification 5 by input 
of I0PF«4. Ihe use of 10F? enables the user to extend the table by adding points 
upstream. 

^ Jf punched cards for input to a boundary layer program are required , the 

fo!lcv.’ing items must be input as pert of the $0D£ Namelist input: 

Item Name Description 

' " IPIAB = If IPTAB*’ , one title card will be punched (this will 

bo the last title card input as described in 6.2} 
followed by tables of X, Y, and P/P along the nozzle 
wall. These cards are for input to BLIMP or 
MABL computer programs. The first point punched will 
correspond to beginning of the converging section of 
the nozzle (i.e. at ECRAT; see Figure 6-1, also 
table 6-8). 

10FF s The first point to be punched will be numbered as 

I(?FF+1. 

■' IPTBL » If IPTBL and IPTAB-1, tables of M, T/T , C , V, and 

R p, will also be punched. These additional cards are 

^ required for input to the TBL computer program (i.e. 

the December 21, 1973 version). 

If a TDE problem is specified, the following items must also b'i input 
' when IPTAftsl: RSI, RWTU, THETAI, and RI in $DATA. 


These items define the nozzle geometry from the combustion chamber through the 
throat as shown in Figure 6-1. For a TDE option it is necessary that IPTBL«o. 



6.5 


OOK INPUT DATA 


0OK input data is required if 0DK =1 or TDK =1 in $DATA as described 
in Section 6.3. The OOK input data consists of three data groups as follows: 


SPECIES 
REACTIONS 
$ 0DK 


data group 
data group 
data group 


These data groups are described below in sections 6.5.1, 6.5.2, and 6.5.3, 
respecti vly. 


6.5.1 SPECIES 


Species used by the computer program are determined In several possible 
ways, depending upon the problem type. Methods used to determine chemical 
species for each probelm type are discussed below. 

m. 

For 0DK problems species names and concentrations must be input, see 
Section 6. 5. 1.1. 

GfDE-GfDK 

For C^DE-Gl’DK problems the initial start conditions for the kinetic exi,':’- 
sion are obtained from an equilibrium calculation. The species list generated by 




the equilibrium calculation generally contains many more species than the 40 
species for which the (21DK subprogram Is dimensioned. Therefore a selection 
processes is required to Interface the 0OE calculated equilibrium start conditions 
with the 0DK kinetic expansion calculations. Ihis selection is performed using 
tiie following rules: 

Rule 1 If a species appears In a reaction. It is selected for the kinetic 
calculation. 

Rule 2 If a species is specified using INERTS directive it is selected 
for the kinetic calculations. 

Rule 3 If any species has a mole fraction greater than an input criterion, 
it is selected for the kinetic calculation. 

Species which are selected but which do not appear in a reaction are treated as i 

Inert and listed as such on the output list of selected species . 

GBDE-GfDK-TDK 

For 0DE-0DK-TDK problems species are selected by the above rules for 
0E>E~0DK problems . However, for multizone TDK cases It Is necessary that 
each zone have the same species list. Ihus the INERTS input (see Section 6. 5. 2. 7) 
must >'<e used to assure the same species are selected for each zone. 

TDK with Input Initial Line Option 

i 

Species names and concentrations must be input as described In ; 

Section 6.8.2 when TDKIL -1 is input in $DATA. 

6.5.1. 1 0DK OPTION FOR INPUT OF INITIAL SPECIES CONCENTRATIONS (APPLIES 
ONLY TO THE 0DK PROBLEM TYPE) 

This input begins with a single card with SPECIES in columns 1 through 7 
and with either MASS FRACTIONS or M0LE FRACTIONS in columns 9 through 22. 

If the identifier for mass or mole fractions is omitted, mass fractions are assumed. 

Up to 40 species cards may be Input. OtJy those species specified by input 
species cards will be considered for an 0DK problem . The order of the Input 




species cards is Independent of the order in which the species appear on the 
master Ihermodynamic Data file. 


A chemical species Is Identified symbolically by 12 alphanumeric char- 
acters and must correspond identically with the species name as it appears 
on the Themoc.mamic Data file. A complete list of the current species names 
are listed in Table 6-4 (condensed species, however, may not be specified in 
the species list.) The species symbol may not contain the characters * or 


Col 

1-10 

11-22 

23-30 

31-60 

61-80 


Function 
Not used 

Species symbol Oeft Justified) 

Not used 

Value of initial species concentration (if zero 
must be input as 0.0) free field F or £ format 

User Identification if desired 


6.S.2 REACTT0NS 

Chemical reactions must be input if 0 dk ■! or TDK »i in $data. 


Up to 50 reactions with an implied third body and a total 150 reactions may 
by input to the pr^^gram. Only one card per reaction, and only one reaction per 
card is permitted. Cards specifying third body reactions must precede cards 
specifying all other reactions. Species names appearing in the symbolic reaction 
set must correspond identically with the species names as they appear in the 
master Thermodynamic Data (see Table 6-4). A card listing for a sample reaction 
set is presented in Table 6-10. 

The symbolic reaction set contains directive cards and reaction/data cards 
in groups as outlined below: 

REACTT0NS Directive for start of symbolic reaction 

input 

« 

, Reactions with implied third body speices 

END TBR REAX Directive for end of third body reactions 

• 

. All other reactions 
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LAST REAX 
INERTS 

THIRD B0DY REAX RATE RATI0S 


Directive for end of reactions 

Specified Inert Species 

Directive for start of third body reaction 
rate ratios 

Third body reaction rate ratios 


LAST CARD Directive for end of REACTIONS input 

The content and format of each type of card is defined as follows: 


6.5.2.1 The symbolic reaction set begins with a card containing the word 
REACTI0NS in columns 1 through 9. Other columns on this card can be used 
for comments . 

6. 5. 2. 2 Each card defining a reaction is divided into five fields, separated by 
commas. Each field contains: 

field I a reaction 

field 2 As followed by a value of A 

field 3 Ns followed by the value of N 

field 4 B = foHov.’sd by the value cf B, 

the activation energy (Kcal/mole) 

field 5 available for comments 

The general fonn of a reaction is: 

Nj^ *Symbolj + N 2 *Symbol 2 ^ * * “ ^a ^ 

where the left hand side represents reactants and the right hand side represents 
products. The reaction can be either endothermic or exothermic. 

The multipliers, N, must be Integers and represent stoichiometric co- 
efficients. If no stoichiometric coefficient Is given, the value 1 is assumed. The 
dimensioning currently in the program requires that: 

N j + N 2 + . . . . s: 10 
and 

N, + N. + . . . . s 10 
a o 

The chemical species (denoted by the word “symbol" in the above general 
form) can contain up to 12 characters each of which must match a species name con- 
tained in the thermodynamic data (see Table 6-2, card 3). 


rate parameters 
for the reactions 
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examples: 

Reaction 

NA++CL-»NACT, 

B+2+0-2«B0 

BE+2+2*0H-»BEC?H0H 


Interpretation 
Na”^ f Cr- NaCi 
B'^ + CT' * BQf 
Be^ + 20H" = Be(C?H)2 


The value assigned to A, N, B define the forward (i.e. left to right) reaction rate, 

k » A • • ^-(lOOOB/ftT) 

In units of cc, ^K, mole, sec. 

All three reaction rate parameters must be input. The numeric value of each para- 
meter may be specified in either I. F, or £ format. If £ format is used the £ must 
appear before the exponent. 


6. 5. 2. 3 The reactions with an implied third body must precede other types of 
reactions, and must be followed by the directive (columns 1 through 12): 

Col 1 ^ 

C?;D TSR R£AX 

all reactions prior to the above directive will have a third body term added to each 
side of the reaction. £.g. 

H2 » K + H, . . . ■ 

£ND TBR R£AX 

defines the chemical reaction 
H2 + M = H + H + M 

where M is a generalized third body. Specific third body effects may be included 
by Inputing specific third body reaction rate ratios as outlined in 6. 5. 2.8. Cards 
encountered after the £ND TBR R£AX directive card do not have a third body term 
added. 


All other reactions are input next, exactly as described under 6. 5. 2. 2. 

6.5.2. 4 After the last reaction has been defined, a card with LAST R£AX in 
columns 1 through 9 is input. 
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6.5.2. 6 Reaction rate data for 13 dissoclation*recomLlnatlon (Implied third 

bod/) reactions and 35 binary exchange reactions are listed in Table 6-10 for 

propulsion systems containing dements C, Cl, T, N and O. These rates are 

% 

taken from Reference 16 ^o additional reactions from Reference 15 are included). 
Cards can be abstracted from Table 6-10 for input to the computer program. For 
the Implied third body reactions, the third body for which the rate applies is in- 
dicated in parenthesis in the comment field (M represents a "generalized" third 
body, see Section 2 2 for further details). 

6. 5. 2. 6 INERT SPECIES OPTION 


Inert species (i.e. species not appearing in reactions) can be included in 
the input by input of a card with INERTS in columns 1 through 6 followed by a list 
of inert species names. The species names must each be followed by a comma 
and each name must written exactly as in the master Thermodynamic Data. The 
last comma must be followed by the word END. See Table 6-11 for an example. 
The species list can continue on to the next card, but a species name can not 
overlap onto Uie next card. 


6.5.2. 7 THIRD BODY REACTION RATE RATIOS 


th 


As described above in Section 6.5. 2.2 for the J reaction only one reaction 


rate, k^, where 


k^ = AT“^J ' 


can be input. For three body recombination reactions such as 
H + C^H + Mj = + Mj 

the rate of reaction is in general different for each species, M^, depending upon 

the efficiency of the species, M^, as a third body collision partner. As discussed 

In Section 2.2 the temperature dependence of a recombination rate is approximately 

th th 

independent of the third body, i.e, or the i third body and J reaction: 

k^rAyT-Nj.-Oj^T 

The third body efficiency of ^he 1^^ species for the reaction is then defined as 
my eAy/A, 
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TABLE 6-10 REACTIONS AND RATE DATA FOR C.Cl.F,H,N. AND O SYSTEMS 
CFROM REFERENCE 16) 


»rAfTrrtw« ccLFHMn h« 

•«»A 197? 

J4NNAF FSWa KaAT*4-N 

EXPI-10A08/RT) 

1 

H ♦ 

r»M 

H70 

••7.5F?3 


Na?.6 t 

«««.• 

(AR) 

MO, 

0 ♦ 

H 

OH 

AaA.OFlA 


w«l» • 


(AR) 

NO. 

2 

0 ♦ 

n ■ 

0? 

<al.2F17 


Mai, , 

nao.« 

(AR) 

MO. 

3 

r ♦ 

f 

P? 

AaS.TElS 


M«l. • 

«*0,» 

(AR) 

MO, 

4 

H ♦ 

P 

HP 

AaE.SElA 


Nai, , 

Rbq,, 

CAR) 

NO. 

9 

H ♦ 

M m 

HP 

A>«.*F17 


Mai, • 

BBO,* 

(AR) 

NO. 

6 

ro? 

0 ♦ CO 

A»2.7F3? 


Ma4,5 • 

B*127,5SS, 

(HP) 

MO. 

7 

r ♦ 

0 

CO 

Aa3*o€t6 


Na,S • 

Rbo, , 

(M) 

MO. 

B 

M » 

N 

H? 

Aal.OEl* 


Nal. « 

RbO.1 

(H) 

NO. 

9 

N ♦ 

n 

HO 

«*6.4F]A 


Ma,;i • 

BBfl,* 

BAULCH (Np) 

MO. 

10 

ct ♦ 

r 

CLF 

Ab3*0E|6 


Mb, 5 • 

Mao,, 

(H) 

NO. 

11 

H ♦ 

a 

HCL 

Aa3«aei« 


Ma,5 * 

Mao, • 

(M) 

NO. 

12 

CL ♦ 

CMO THIS 

CL 

} OPA 

CL? 

A»i.iei* 


Mai, • 

B«0*« 

(M) 

NO. 

13 

H? 4 

OK 

H * H?0 

A«2.i«ei3 


MaO, 1 


BAULCH 

NO. 

21 

OH ♦ 

OH 

0 * H20 

A>5.7sei2 


Mao, • 

Bb,7Bo* 

RAULCH 

NO. 

22 

H ♦ 

OH 

0 * M2 

Aa7.33Et? 


Mao, « 

Rb7,3oO* 

BAULCH 

NO. 

23 

0 ♦ 

OH 

H ♦ 02 

Aal,3Fl3 


Mao, • 

B««.» 

BAULCH 

NO. 

25 

OH • 

cn 

H ♦ cop 

AaS.Afll 


Mao, » 

Sai.OROt 

BAULCH 

NO. 

Al 

02 ♦ 

CO 

0 ♦ CO? 

Aa«.ASC9 


Ma.,ASA,Ra49,9?0, 

BAULCH 

NO. 

A2 

H ♦ 

NO 

0 ♦,H2 

«a3.iFi3 


N«;i» 9 

BB.3i4, 

BAULCH 

NO. 

4B 

H ♦ 

OP 

0 * HO 

Aa6.43C9 


H«»l« 9 

9>6,290* 

BAULCH 

NO. 

51 

OH « 

OH 

M2* 02 

* *al.4)ei3t 

Ns«o159 

Ba49.?64, 


NO. 

24 

H ♦ 

P? 

HP* P 

t 4aS«3F]2 

9 

N»»,S 9 

9*A.0e0* 


NO. 

25 

H2 ♦ 

P 

HP* K 

* taS.OEl? 

« 

Hv^« 9 

BaS.700* 


NO. 

27 

M7 ♦ 

P? 

HP* KP 

t Aal.7R£lO( 

H«-,5 9 

8«39.73»» 


NO. 

2B 

H ♦ 

a? 

HCL*CL 

» *a3.()Fl4 

• 

9 

B>3,00A« 


NO. 

29 

♦ 

a? 

HCL*HCL 

• 4al,7A£lo» 

N«»*S 9 

9a45.37S 


NO. 

30 

HCL* 

H 

Hp*CL 

• *s6.2Ft1 

9 

N*»,5 9 

B«3.{00* 


NO. 

31 

MCL* 

P 

HP*CL 

. Aal,9Fl2 

• 

Hs*«6S9 

BMfBnA* 


NO. 

32 

CL?* 

P 

CL*CLP 

t 4a«,?Et? 

9 

Nh-9699 

Bb,5AA* 


NO. 

33 

CL ♦ 

P2 

p *cLr 

, «a7.6Ej2 

• 

Hs»«f09 

BM*30Af 


NO. 

34 

CLr* 

H 

HP*CL 

* Aai.8Fl2 

9 

He«»609 

9M3«2A0t 


NO. 

35 

CLr* 

H 

HCL*? 

• 4aS.6Cl2 

9 

N««m«M9 

0m1,9oo« 


NO. 

35 

cir* 

H? 

hCL*hp 

» Aai.SFlo 

t 

9 

0b4A,337, 


NO. 

37 

r? ♦ 

HCl 

HP*CLP 

• AaltAClO 

• 

N«*«9 9 

0a3B.*27, 


NO. 

3B 

CL^* 

HCL 

HP*CL2 

• 4al.8EtO 

9 

Nb»«9 9 

0b*6,o2*» 


NO. 

39 

f? ♦ 

CL? 

CL?*CL? 

t Aal.RElo 

• 

N«*.S 9 

8B26.790, 


NO. 

AO 

CO?* 

C 

co*co 

• F»l»lEll 

• 

H«-t5 9 

0a6.9Q9, 


NO. 

43 

c ♦ 

OH 

CO* H 

• 4S5.3F11 

9 

9 

0a9.A7B, 


NO. 

44 

C ♦ 

NO 

CO* N 

t iaS.SFli 

9 

N«-,5 9 

0bB.3o3* 


NO. 

49 

Cft?* 

N 

CO*HO 

* Aal.lEll 

9 

HH»«9 9 

0aS4.«10, 


NO. 

46 

c ♦ 

0? 

CO* 0 

• 4a5,3En 

9 

M»-.5 9 

SaO.SB?, 


NO. 

47 

HO ♦ 

NO 

HP*02 

« *al.0F13 

9 

H»S. 9 

0a79,4OO, 


NO. 

49 

H ♦ 

OH 

HO* H 

, 4aS.3F11 

9 

9 

0aS.A?A, 


NO. 

5o 

OH ♦ 

P 

HP* 0 

. 4B2.9E12 

9 

H«»«609 

OB*?OAt 


NO. 

92 

H20* 
HCL * 
OH ♦ 

last pi 

P 

OH 

CL 

FA< 

HF*OM 

H20*CL 

HCL*0 

• «al.4Flo 

• Aai.nFn 

• taS.OCll 

9 

9 

9 

Nv»»909 
H«*,9 9 

nb«,9 9 

0a«OoA, 

0B«.O* 

0mA,O* 


NO. 

93 
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Thus m^j is the ratio of the reaction rate with species as the third body to 
the reaction rate input on the reaction card described in Section 6.S.2.2. 

If reaction rate ratios, m^j, are to be input for the dissociation-recombin- 
ation reactions, a card with THIRD B0DY R£AX ItATE RAUC^S in columns 1 through 
27 must be input next. If this card is deleted from the input, the program assumes 
all m ^^^ » 1 . If this card is included in the input, it must be followed either by 
a card with ALL EQUAL 1.0 in columns 1 through 13 (which sets all ot by 

SPECIES cards as described below: 

The m,(, can be input using a card with the word SPECIES in columns 1 through 
U th 

7, Ihls word is followed by the name of the 1 species followed by a comma, 
followed by the values mj^j in F format, each followed by a comma. These m^j values 
can be continued onto succeeding cards . Note that the m^^ values depend on the 
order of input of the reaction cards, i.e. the reaction is defined by the 
card input after the REACTIONS card. 

Table 6-11 gives a sample input for a Hydrogen/Oxygen system using third 
bpdy reaction rate ratios. In this example the three body recombination rates 
are input with Argon as the third body. The rate with H 2 as a third body is a 
factor of 5 larger then with Ar as a third body for the first three reactions and a 
factor of 4 larger for the fourth (Hydrogen recombination) reaction. 

6. 5.2. 9 At this point in the data input deck a card with LAST CARD in columns 
1 through 9 must be input. 
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TABLE 6-11 USTING OF SAMPLE REACTIONS CARDS FOR AN 0,/H, PROPELLANT 

4b (6 


REACTIONS 

D-H 

MAY 

3-4 1072 , 

lAi^ikiAr 

PSWS 



H ♦ OH 

« HEO 


• 

A«7*5E23 • 

N*2,6 

• B>0«t 

(AR) 

NO. 

0 ♦ H 

■ OH 


• 

A«4*0El8 • 

N«l, 

. t >6 bo« • 

(AR) 

NO. 

0*0 

• 02 


• 

A-U2E17 • 

N«l, 

t B«0»t 

(AR) 

NO. 

H • H s H2 
END TBR REAX 


• 

A«6,4El7 , 

N»l, 

• B*Ott 

(AR) 

NO. 

H? « OH 

• H ♦ 

H20 

'/ 

A«2.19Ei3t 

NaO, 

• 8«9.1St 

BAULCH 

NO. 

OH ♦ OH 

■ 0 ♦ 

H20 

f 

A«5«75Ei2* 

H»0* 

• 8«.7Bot 

BAULCH 

NO. 

H ♦ OH 

■ 0 ♦ 

H2 

f 

A«T,33E12# 

N*0* 

• 8*7»3oO* 

BAULCH 

NO. 

0 ♦ OH 
LAST REAX 

• H ♦ 

02 

• 

A«lt3El3 • 

N«Ot 

t 8*0 •• 

BAULCH 

NO. 


{NmS N9«AA»EN0 

THIRD ROOT R£AX RATE RATIOS 

SPECIES AR«U«U»l.tl*» 

SPECIES H?tS«tS««S,f4.» 
species NaOt20,*3«,s«*?6.« 

SPECIES 02*S.*S,»4.StUS« 
SPECIES N2«4,$4.t4,«U5« 
SPECIES H,l:>,Stl2«5«)2.S«2S,t 
SPECIES 0tl?*5«12«5*l?*5«2S»t 
SPECIES 0 H«i2,S,i2«s»i2«5«25,» 
LAST CARO 



6.5.3 $0DK NAMELIST INPUT . 


$0DK Namelist Input specifies the conditions for the kinetic expansion calcula- 
tion. The input is read in subroutine ^'DKINP and consists of the following groups 


SMJLSJfJk CAk» WC* IaXJUA tfCVA * 


6.5»3*1 Reification of Nozzle Geometry 

6 . 5 . 3*2 Integraticai Control 

6 . 5 . 3*3 Print Control 

6 . 5 . 3«4 Rcies Selectiai and Mass/Mole Fraction 

Qieck 

6 . 5 . 3*5 0DK Problem Input 


6 . 5 . 3.1 Specification of Nozzle Geometry . 

All of the nozzle geometry is to be Irput using the $DATA Namelist irput, see 
Section 6.3.3 and Figure 6-1. 


An JWK calculation will be carried out for each mixture ratio ir^ut in 9 FSKED 
of $0DE whenever (2)DK = 1 or TDK » 1 in $DATA. 

For a TDK problem, it 1s necessary that tne 0DK calculations be run past the 
nozzle throat. Usually it is not desirable to run the 0DK calculations all the way 
to the nozzle exit because the extra computer time and print out that results. 
However, if this is desit-ed on a TDK problem, it can be requested by input of item 
EP as described below. 


Item 


Description 


Assumed 
Units Value(s) 


$0DK 

-s 

Namelist, read in subroutine 

0DKINP. 

EP 

3 

If TDK » 1 and a value is input 
here for EP, then each ODK will be 
run to expansion ratio EP. 
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6. 5. 3. 2 INTEGRATION CONTROL 


The intenraticr. routine controls tho step size such that the relative 
error in the dependent variable increments are less than a prescribed fraction, 
DEL. Only doubling or halving of the step size is permitted, and on option, 
either all the variables may be considered (JF=0) , or only the fluid dynamic 
variables (JF=1) may be considered. 

When the flow becomes supersonic and the area defined fluid dynamic 
equations are used, an additional check, on continuity is applied in the form 


(pVA)^^^ - (PVA)^ 


< C<2^NDEL 


where C0NDEL is an input relative criterion. 

The step size is held between the two input bounds HMIN and HMA.X. 
Fixed step cases may be run by setting input values for KI, HMAX, HMIN all 


equal . 
Item 


Description 

UnU^ 

Assumed 

Value(s) 

$0DK 


namelist, read in subroutine 
ODKINP . 



HI 


initial step size 

non^ 

.01 

HMAX 


upper bound on step size 

none 

0.10001 

HMIN 


lower bound on step size 

none 

.005 

DEL 


fractional incremental 

none 

.001 



error 





V 

1 
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Item 

TEXPLi 


C0NDEL 


IF 


0 


Description Units 

temperature below which 
explicit integration will 
start . Not recommenrlpd . 

relative error criterion for none 

continuity check for super- 
sonic flow 

all variables considered none 

for step size control 

only fluid dynamic variables none 

considered for step size 
control , i . e . , p , u , and T 


Assumed 
Value(s ) 

0 


1x10“^ 


0 


i 






6. 5. 3. 3 PRINT CONTROL 


Output from .the Kinetic Expansion Calculation consists of complete 
output for each print station selected. The end pomt of the nozzle is "Iways 
piinted. Piint .«tatidns are selected from one of the following options: 

Assumed 

Item Description 


JPRNT 


-2 print throat and input area ratios (see ARPRNT) -2 

-1 print et internally set a>~ed rati os for conical 

nozzle.* Print at selected wall contcur pointc 
for contoured riozzles . For t‘ t sp.ine fit- 

option (IWALL=4), print out will occur at 
each ent^y in ZS of $3PP. For other con- 
tours (IWALL = 2 ->r 3 )* pri^c out will 
occur at 20 equall;> spaced axial locations 
along the nozzle. 

0 print at every integration step 

+1 print every ND3rd step up to the throat and 
then nominal area ratios 

+2 print every ND3rd step over entire nozzle 


I • 



*For JPRNT ---1 a.id a conical .nozzle (1 FA^ALL = 1), the internally set a.rea 
ratios arc: 

A.RPRNT(l) --- 2,3, 1, . . 39,40,42, , . , 58 , 60 . 64 , . . . 1 16 , .1 20 . 128 ...,200,210, 

220,. .400 





i 
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Item 


D escription 


NDi 

ND2 

ND3 


first integration step to be selected for 
print 

last integration step to be selected for 
print 

print every ND3rd step between NDI and 
ND2. 


If JPRNT is -2, the following must be input: 


Item 


Description 

ARPRNT(l) 


requested area ratios for print, must be 
monotonic increasing and greater than 1.0 
If no values are input, will use 
values from ASUP of $SPP. 

N JPRNT 

=: 

number of area ratios requested for orint 
< 100. 

extended print option may be selected as follows: 

Item 

Value 

Description 

IDYSCI 

= 0 

no extended print requested 

• 

= 1 

extended print option selectee (not suggested) 



OR!QW'!AL 
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6. 5. 3. 4 SPECIES SELECHON AND MOLE/l^SS FRACTION CHECK 

In order to Interface 0DB equilibrlunfi calculated start conditions 
with the kinetic expansion calculations, special consideration must be made 
for inert species (those not appearing in the reaction set). Inerts may be select - 
ed explicitly by use of the INERTS directive or by use of a relative selection 
criterion. 


If A MULTIZONE TDK PROBLEM IS SPECIFIED INERTS MUST BE SPECIFIED 
VIA THE INERTS DIRECTIVE. This is required so that the chemistry selected for 
multizone cases will be compatible. 

The INERTS directive is described in Section (6.5. 2.6). 

The relative selection criterion (0DK or 1 Zone TDK prdtjlems,) 
is described below: 

Item Name Function 

EPSEL * sll species which do nor appear explicitly 

in the reaction set but whose mole fractions 
are greater then the input value for EPSEL, will 
be retained for the kinetic expansion. Species 
selected under this criterion are treated as inert. 
The program assurries £PSEL = l.OE-5, unless input. 


In some instances it may be desirable to use input species concentrations 
which do not sum to unity. Species concentrations, either input or from equili- 
brium start conditions, are summed and the sum checked as described below. 


Item Name 


Function 


XMFTST ** Input species concentrations are summed and checked 

versus unity using this input crite:ion. If 

j 1 - J^species concentrations j<XMFTST 

then the test is passed. The species concentrations 
will then be noTr.alized such that 
2])species concentrations = 1. 

The program assumes XMFTST = i.OE-3, unless input. 

If the test is not passed, an error message will be given and the 
run terminated. 
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6.5. 3.5 ODK PROBLEM INPUT 

This input is required when PROBLEM 0DK is specified on the 
problem card. A kinetic expansion from input arbitxdiy start conditions is 
to be computed. In addition to the input items described in section 6.5.3, an 
0DKproblem requires input of those items described in sections 6.5.1 and 2. 


Item Name 


In out Quantity 

Units.. 

SI Units 

PC 

- 

chamber pressure 

PSIA 

N/M* 

T 

» 

Initial temperature 


®K 

V 

- 

Initial gas velocity 

ft/sec 

m/sec 

JPTIAG 

m 0 

pressure table 
calculated internally 

none 

none 


- 1 

pressure table input 


« 

ECRAT 

m 

initial conti action ratio 

none 

none 

For JPFLAG * 0 option the 

following must be input: 



Item Name 


Inout Quantity 

Units 

SI Units 

PI 

a 

initial pressure 

PSIA 

NAI* 

PESTAR 


throat pressure 

PSIA 

N/M* 

For JPFLAG * 1 option the followii, 

■g must be input: 



Item Name 


Incut Quantity 

l^nits. 


PTB(l) 

sr 

normalized pressure 
table entries* 

none 


ZTB(l) 


normalized pressure 
table coordinates** 

none 


NTB 

“ 

number of pressure 
table entries, s 127 

none 

t 

Z 


initial axial position 

none 



* normalized to Input chamber press'ore, PC 
** normalized to input throat radius , RSTAR 
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6. 5. 3. 6 MASS AVERAGED g)DK BP 

A mass averaged (&DK ISP summary page may be obtained at the end of 
the 0DK calculations as described below: 

ftgJP Ngmg Pescrlpfioo 

MAVISP *1 Specifies mass averaged ISP option 

XM(1) s Ratio of mass flow rate of each zone (zone 1 

* inner zone) to the total mass flow rate. 



6.6 


$TOANS NAMELIST IIPrJT. 


Wien a f<8C problem has been specified, the input data set $TOANS is required 
for the transonic calculation. 


Item 


Description 


Assumed 

Value(s) 


$TRAN 

3 

Namelist, read in subroutine TRAN 

50*0 

XM(1) 

S 

Ratio of mass flow rate of each zone (zone 1 
» inner zone) to tiie total mass flow rate, 
(need not be irput if MAVISP = 1 cation spe- 
cified and XM ir^t in $0OK). 


ALI 

3 

lAttDber of degrees initial line will be disp- 
laced downstream. Ihe program assumes ALI 
is zero. If ALI is not zero, a symmetric 
throat is required (RWTD » RWTtJ). 

0 

IBUG 

3 

If input is nonzero, intermediate transonic 
output will be printed. 

0 


Tlie following input way be used to control the construction of the initial line: 


Assumed 

Item Description Value(s) 


MP Number of points to be placed on the initial 50 

line. For RWTD <_ PMdLT the number of 
points placed on Bie Initial line will be 
Incremented by the largest integer of theta, 
the downstream tangent angle. 

PMCRIT For RWTD > FMCRIT a sinusiodal distribution 1. 

is requested; i.e. 


i TT 

r^_ « [r^sin (|j 2 ^] i = 0,1,2 N 

wtjre N = MP and e-j is EXP1 descrioed below - 
Editing is dcMie to control the spacing. 

If RWTD <_PMCRIT an initial line distribu- 
tion will Fesult such that the first left 
running characteristic intersects the wall 
at a wall angle of PMDEXjr degrees. 
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Item 


Description 

Assumed 

Value(s) 

PMDB3 

a 

Intersection angle, in degrees, for tiie 
first left running r.haracteristic if R^TD < 
PMCRIT. 

1. 

EXP1 

3 

for sinusoidal distribution. 

1.2 


$END 



6.7 


$MOC NAMELIST INPUT 


This data set contains the input items for the supersonic Method of Character- 
istics (MOC) module. The items are divided into four types, which are described 
in the following subsections. 

6.7.1 Characteristics Mesh Control 

6.7.2 Inputs from SDE3R 

6.7.3 Print Control 

6.7.4 Exit Plane Option 

Often no $MCX^ input is necessary since the default values are usually 
ficient. 


suf- 



0 


^ 

■r 


'{ 

:| 1 

1 6.7.1 Characteristics Mesh Control 


1 

Item Name 

Input Quantity 

Uni ts 

j 

tMnr - 

IV V 

Namenst name, read by subroutine 
CHAR 


! 

DS 

Insertions will be made such 
tiiat successive points along 
streamlines will not be sep- 
arated by more than DS. 

none 


DWWI 

Insertion control parameter 

described In Section 5,9,2, 

degrees 


EPW « 

the program will insert such 
that the wall end point is 
located within a tolerance EPW. 

none 


IMAX 

the maximum number of interatlons 
to be allowed while attempting to 
achieve a relative convergence for 
the flow variables of 5*10“S. 

none 

r 

IMAXF « 1 
or 

IMAXF - 0 

tiie program will terminate the case 
if a printed point requires maximum 
iterations for convergence. 

program will continue the case after 

none 


IMAX iterations per pclnt have occurred 

none 

1 

TEXPU » 

input temperature below which explicit 
integration for the species concentra- 
tions will be used. 

*K if SI Units 

• 

ETHI 

(gfor point editing as described in 
Section 5.9.2, CNTRL. 

degrees 

1 

ES 

t ^ for point editing as described in 
Section 5.9.2, CNIH.. 

none 


DIWI * 

for insertion in sub- 
routines INPT, DSPT as described 
in Section 5.9,2, CNTRL. 

degrees 


i 

1 
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6.7.2 Print Control 




i 

i 


1 


.1 

i 

I 


Itsrn NomS 
N1 


N2 

NC 


h^ASSFL 


NDS 


input Quantity 

flow parameters will be printed 
for every Nl^ Interior point a- 
long characteristics selected 
for print 

every N2^ characteristic will 
be selected for print 

for NC ^ 0 species concentra- 
tions (partial densities) will be 
printed with the flow parameters. 

If NCsl the quantities 
heat capacity (BTU/Ib-®R), and 
enthalpy (ft*/sec*) will be append- 
ed to the species concentration 
print. 


at the completion of each left running 
characteristic (LRC) the massflow is 
integrated. 

If 

MASSFL » 0 then no mass flow printed 

MASSFL = 1 then total mass flow and the 
number of points on the LRC 
are printed for each LRC 

MASSFL « 2 then mass flow for each point 
along LRC is printed 


MASSFL « 3 Same as MASSFL = 2 with the 
addition of execution time at 
the end cf each LRC 


see Section 5.8.1, CHAR 

for NDS = 1 Dividing Streamline Points 
will be printed. (Nominal) 

for l^S s 0 Dividing Streamline Points 
will be suppressed. 


Units 

none 


none 

none 


none 


none 


I 
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6.7.3 Inputs from DER, Reference 10 




Item Name 
0FBAR 

ETABAR 

DRPISP 


Units 


Overall mixture ratio in- none 

eluding condensed phases. 

For print out only. 


Overall evaporation efficiency, 
l.e. the ratio of gas flow to 
total propellant flow at the 


throat, I , total 
♦ETABAR P 


^sp,gas 


none 


Ratio of total condensed Ibf sec/lbm 

phase momentum to the mass (if SI Units then 
flow at the throat. Not used, N secAg) 
reserved for future use. 




V 


6.7.4 Exit Plane Option 

On option, the TDK method of characteristics calculation will continue 

the mesh construction through the exit plane of the nozzle and print a summary 

» 

of the exit plane properties. 

Item Name Input Quantity Function 

EXITPL ■ .TRUE. Exit plane calculation 

requested. 


6.7.5 Punch Initial Line 


During any calculation generating an initial line, the initial line 
may be punched in a form suitable for running an input initial line option. The 
following input Is required. 


Input Quantity Function 


UPCH * .TRUE. Requests Punchlnj of 

Initial Line 

IPUNIT The Fortran unit number 

assigned to the PUNCH, 
(e.g. on Univac 1108 
IPUNIT « -a, on the 
other machines it may be 
7, 8, etc.) 


IMPORTANT NOTE 

IF IPIAB » 1 option is selected, i. e. the boundary layer edge conditions 
punched f'^r TBL input, the initial line punched cards will be Interspersed 
with the TBL edge conditions punched cards. 


SEND 



6.8 


Boundary Layer Module (BLM) Input Data, $BLM 




This input data set is required if the BLM is to be executed, i.e., if BLW=1.0 
was irput in the $DAXA namelist. Most of the data required by the BLM is cammijnica- 
teiu auwjmafcicaiiy by the 0CE and/or M0C modules, or is preset as assumed values. 
Hence, these data itans do not necessarily need xo be irqjut to the BLM module. 
However, any value that is read in will override the assumed or communicated value. 

The ir^jut data items to the BLM modulo are as follows. 



Item 


Description 

t&iits 

Assumed 

Value(s) 


$BLM 


NEunelist name, read in Subroutine 
IMPUTE 



r 

ITYPE 

= 

flag to specify the type of body 
geometry 

1. 


i 



for an axisymmetric nozzle, re- 
quired by SPP 



r 

> 

f 


- 2, 

for an axisymmetric external flow 



.* 

c 

\ 


= 3, 

for a two-dimension?.l external 
flow starting at a stagnation 
point. 



i, 

WD0T2D 


nozzle mass flov’ If the TD2P 

Ib/sec 

0. 


module was not e> .ecu ted, a value 
can be irput here so that a boun- 
dary layer ISP decrement can be 
computed 
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6 . 8.1 


ELM Gas Properties 


Gas properties that are required by the BLi'l ar^ y, Cp, and u versus T, and a 

value for the Prandtl number, P„. 

r 

If tables for Cp and Y are not input (.‘Joe CPO, CK), and TO, below) then the 
program will prepare these tables usLng the 0M) module. The tables are prepared 
using a series of (T,S) equilibrium calculaitriV., where T varies from 600°R to 
7000°R at 20(5^R increments. The chamber j-y tr-of-v value is used for S. Values at 
100°R are then extrapolated and added to the ta>'ie. The !>able is printed with the 
BLM output. 

If constants defining the gas viscosity, u , and if the Prandtl number, P^., are 
not irput, then they will be transmitted from the ODE module. 

If the TD2P and 0DE modules have not been run, then the data listed below must 
be irput. 


Item 


Description 

Units 

Assumed 

Value(s) 

CK0(1) 

s 

table of ratio of specific heats, 
Y, versus T. 

none 

101*1.4 

CP0(1) 

3 

table of specific hrat at constaiit 

pressure, C_, versus T 
F 

ft/sec-R 

101*6006 

TO(1) 


temperatures corresponding to cbe 
entries in CKD and CPO arrays, 
above 

^R 

300,. ..,7000 
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Item 


Description 

Units 

Assumed 

Value(s) 

•“ 

NTAB 


nuiJiber ox values entered in CKO, 
the CPO, and TO arrays. 3 < 101. 
< MTAB 

- 

24 

'j 

RMUI 

sc 

Referep'^ viscosity, y^, where vis- 
cosity is expressed as 

Ibm/ft-sec 

.605934*10”^ 

1 

1 

} 

<1 






1 

TI 

* 

reference ten^jerature, T^, for vis- 
cosity. See RMUI, above. 


6032 


0MEXIA 

3 

Viscosity exponent, 0 ). See RMUI, 
above. 

" 

.66843 


PR 

3 

Molecular Prandtl nuniier, 

- 

.4385 




6.8.2 Boundary Layer Edge Conditions 


Ihe coordinates for the boundary layer are specified in the RIN0 versus XII^ 
table. Conditions at the invisid edge are specifed in the UE50, TE0, and PE0 versus 
XIN0 tables. Conditions at ISie wall are specified in the TQV and CQV versus XDQW 
tables. 'Hie program will redistribute the irput stations (up to 201 total) in 
order to have 101 x-staticHis uniformly distributed per segment, except for the 
fj.rst five stations vrtiich are generated non-uniformly. llie input values of 
(or a ), and (pV)„ are interpolated at the new x-stations and used in 
the boundary-layer calculations. 

Values of RINO versus XINO and the conditions at the invisid edge of the boun- 
dary layer will be automatically transmitted from the TD2P module if it has been 
run. Otherwise, th^ must be input here. See the description of I0FF in $DA1A. 
Condi ticHX at the wall must always be iiput here. 


Assumed 


Item 


Description 

Units 

Value(s) 

XH«5(1) 

3 

x^ axial coordinate 
e 

none 

- 

RIN0(1> 

S 

y^, radial coordinate 

none 

- 

UB0(1) 

3 

Ug, gas velocity at the invisid 
edge 

ft/seo 

- 

TE2)(1) 

3 

T , gas static temperature at the 

W 

invisid edge 

0- 

- 

PB0(1) 

3 

Pg, gas pressure at the invisid 
edge 

psl 

- 

NHN0 

3 

number of items jji the XII^, RIN0, 
UB0, TB0, and PE0 tables. NXIN < 




201 . 



Item 

XTQW(I) 

IHFLAG 


TQW(1) 


CQW(1) 

NTQW 


DeacripUon 

X , axial coordinate for the TQW 
and (XM tables. 

Flag specifying wall boundary con- 
dition iiput thro»igh TQW array. 

IHFLAG = 0, for temperature 

IHFLAG « i , for heat flux 

T^, wall temperature or wall 

heat flux, depending in IHFLAG. 

For an adiabatic iirall, set IHFLAG = 
1 and all TQW(1) 0. 

To convert (BTU/hr-ft^) to 

(Ibm/sec^) multiply by 6.9547, 
\diich is Jg/3600. 

(pV)^, mass transfer parameter at 
the wall 

Number of axxal stations, X^. Each 

of the above tables must have this 
number of entries. tfJXM <101. 


Units 


o 


R 


or 3 
Ibm/sec 


Ibm/ft^-sec 


Assumed 
Value (s) 


1 ' 


101*0 
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6.8.3 Integration Step Size Ctyitrol . 




Description 

Unite 

y cr*x uts 

NSEGS 


Itober of Segments, 1_< NSEXJS < 10 

none 

1 

NISPSd) 

a 

Number of Integration St^s per 
Segment, £ 101 per segment. 

none 

10*101 

XSBGd) 

3 

Vector containing the axiad (x) lo- 

ft 

* 


cations vrtiich define the wall seg- 
ments. The vector is always NSBGS 
+ 1 values long and is defaulted to 
the starting and ending values of 
the wall table supplied by the TD2P 
luodule. (The variable XITEL is 
used to define the starting axial 
location for the boundary layer 
edge table in the namelist. 

13ie ending points for the edge tab- 
le is taken as the nozzle lip.) 

NTR s Statical at vrtiich transition to none 3 

turbulent flow is allovred. The 
program starts with an assummed 
boundary layer profile, and then 
turns on the eddy viscosity terms 
for turbulent flow at station ^f^R. 

For a laminar boundary layer, set 
Mill large, e.g. NTR > E NISPS(I). 
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6.8.4 Regererative Cooling Heat Transfer . 


When an engine is cooled using a regenerative device, propellant (usually fuel; 
is routed around the nozzle so that heat is transmitted from the boundary layer to 
the coolant. TMs heat is then returned to the combustion chamber in the form of 
increased propellant enthalpy. When the nozzle wall temperatures are assumed known 
(IHFLA.G a 0 option), the BLM will calculate the heat flux from the boundary layer 

to the wall, fi^(BTU/ft -sec). !hese values can also be input directly (IHFLAG = 1 
option). If the coolant circuit extends from position x^ to position x^ and opera- 
tes at an efficiency, n , then the propellant enthalpy entering the chamber will be 


In) ft (JA (BtU/lUn) 


is the total engine mass flow rate 

is the nozzle surface area differential, 2irrdx 

Using the method outlined above, the BUI will compute increments of propellant 
enthalpy for up to 3 fbel or oxidizer circuits and print cut the resultant enthalpy 
increments. These can be added to a later computer run by using the EELH1(1) iiput 
array. If BU! is to be automatically rerun, then the enthalpy increments will 
automatically be stored into DELHI (1) for the second pass througli OCE, ODK, and 
TDK. These enthalpy increments can be calculated in two ways. If the enthalpy in- 
crease is distributed equally throughout the chamber, then 

= AH 


increased by 


ah ^ 


where 

dA 


l.e. 


DELHI (1) » AH, AH, etc. 
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Ihe second method is to assume that a fuel circuit adds enthalpy only to fuel, 
and an ox circuit adds enthalpy only to ox. It follows that 

for a fuel circuit 


r+1 

4Hi . (^) 4H 

and for an ox circuit 
p*1 

■ <7pr> r 

where r^^ is the mixture ratio of zone i and r is the overall chamber mixture ratio. 

For either method, the steady state engine cycle balance can be approximated as 
follows. First, calculate "adjusted tank enthalpies" for the fuel and for the oxi- 
dizer and irput these on the reactant cards. These values must approximate the 
energy content of the propellant entering the main combustion chamber accounting 
for all energy gains and losses, except heat returned to the main combustion cham- 
ber by the regen cooling circuit(s). Estimates for these amounts are to be altered 
using the E£LH1(1) input array. An estimate of zero is usually satisfactory. Cor- 
rected estimates will be calculated by BIH and stored in IHH1(1) for a second pass 
through TDK (or TI®). A second pass using these values will be executed automati- 
cally if IRPEAT a 1 was itput in the $DATA namelist. 



Item 


Assumed 

Value(s) 




Description 


X0(1) » the Ith entry is the starting position for 3*0 

the Ith cooling circuit*. 

XE(1) a the Ith entry is the aiding position for 3*0 

the Ith cooling circuit*. 

ETA(1) a the Ith entry is the efficiency for the Ith 3*1 

cooling circuit. 

0PC(1) = Type of coolant for the cooling circuit: 3*0 

0FX^(I) =0. if there is no ith circuit 
i®x:;(I) « 1. if the Ith circuit is oxidizer 
^x;(I) a 2. if the Ith circuit is fuel 
I <3 

DISTRBd) a Flag for method of distributing AH incre- 3*0 

malts 

DliraBd) a 0. for equal distribution of 
heat to chamber 

DISTRB(I) a 1. for distribution of ox heat 
to ox in chamber, and fuel 
heat to fuel in chamber. 


■nie Ith entry is for the Ith circuit. \ 


* Normalized by the throat radius, r^.. 





6.8.5 BLM Plotted Output . 


Tlie input described below is used to control plotted output from the BUi. The 
plot options available from BLM are: 

1) Momentum thickness, 9, vs. axial position. 

2) Displacemait thickness, <$*, vs. axial position. 

3) Wall Temperature, T^, vs. axial position. 

4) Velocity profiles at specified area ratios, or at specif isd axial 
locations. 

5) Ten^Derature profiles at specified area ratios, or at specified axial 
locations. 


Example plots are presented in Figures 6»2 through 6-5 for types 1 throu^di 5» 
above. 


Item 

IPR 18 F 

APROF(I) 


NPROF 

KDaPLT 

KMTPLT 

KIVTPLT 


Description 

If “ 0, then supersonic area 

ratios are irput in APRj®'. 

If IPR0F a 1, then axial loca- 
tions, x/r^, are input in APR;®’ 

Area ratios (or axial locations, 

see IPR^®) at which velocity ratio 
and temperature ratio profiles will 
be plotted. Two frames per area 
ratio will be plotted: U/U vs. 

edge 

y/y and T/T vs. y/y 

edge edge edge 

I^ber of area ratios (or axial 

locations) requested in APROF. 
NPROF < 20. 

If iODBPLT * 1, then displacement 
thickness, 6*, vs. axial location, 
X, will be plotted. 

If KMIPLT * 1, then momentum thick- 
ness, 6, vs. axial location, x, 
will be plotted. 

If KTWPLT a 1, then wall tempera- 
ture, T^, vs. axial location will 

be plotted. 


Assumed 

Value(s) 


0 


20*0 


0 


0 


0 


0 


SEND 
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6.9 


SPSCIAL OPTIONS. 


6 . 9.1 C onstant Properties (^s Ont.lnn . 

Ihe TDK program contains a useful option by which the real gas cnemistry can be 
replaced by constant properties chemistry. Multiple zones can be calculated. The 
output includes the nozzle divergence efficic.icy, (see Section 5» subroutine 

PRINT). 

The caistant properties gas option is run by input of PPC30PT 3 1 in $DATA. 
Sample input data for the constant properties gas option is listed in T^ble 6-12. 

Chly the $DATA, $IRAN and $M^ data sets are to be input. These data sets are 
irput as described in Sections 6.3» 6,6 and 6.7» respectively, with the following 
requiroo additions to the STRANS Namelist. 


G(1) * Value of specific heat ratio, Y , for _ach zone, inner to outer, tht 

number of zones is specified in $DATA. 

PSA * Chamber pressure in Ibs/in^. (N/m^ if SI units) 

XPM) = (From Thble 6-9) All assumed = 1, if not irput. 

TC(1) = Chamber temperature, for each zone, inner to outer. ix SI 

units) 

RGC(I) 3 Real gas constant, f+^/sec^ °R, (i.e., 49721 /M^^ 3 g*J*1.986/M^) for 

each zone, inner to outer, (m^/sec^ °K If SI unts) 

XMW(1) 3 Gas Molecular Weight. If input then RGC(1) need not be irput. 
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Table 6-12: Sample Case for the Constant Properties Gas Option. 


TITLE Sample CASE O^C 
DATA 
S0AT4 
FF 60 PT= 1 , 

AIZOmESsI, 

PSIsZ, 

RWTUaZ, R^TOs.5, 

TMETA»35,<»73fl# 

IWALL»4, 

NM 3s 1 1 

R3(2)s* 1 , lo64 3, 1 ,2b«75# I o# i , 7357 S , 2 , OaPaO, <i, 45 ‘» 30 , 

3,b8a26,4.Aa7 72,5,7<Jl98,6,32aoi, 

Z3(2)s ,3<»S75, ,53008, ,8,' ?05, i , 1 9a 73, 1 .66923, 2,32795, 

4.68717,7,68599,10,9601,13.3114, 

THEsU,58l3# 

SEND 

jTRAfJS 

6S1.23, 

PSAalOO, 

TC*5500, 

XM«*20, 

XM»l, 

ALlsO, 

SEND 

IMOC 

"lENO 
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6.10 


INITIAL VALUES FOR THE $0DK. $mNS. AND $TDK INPUTS 

Ihe following defines nominal values to which variables will be set If 
not Input. If a variable Is net-4>lsted, no nominal value is set. Variables are set 
in the subroutine containing the Namelist read.. 

$0DEt set In subroutine 0DES 


DELH(I) 

zs 

0. 

ECRAT 

s 

0. 

EQL 

« 

.TRUE. 

eqthst 

S5 

.FALSE. 

ERATI0 

ss 

.FALSE 

FA 

- 

.FALSE 

FPCT 

» 

.FALSE. 

FR0Z 

8 

.TRUE. 

I{9NS 

= 

.FALSE 

KASE 

= 

0 

LISTSP 

« 

.FALSE 


8 

.FALSE 

CJTFSKEDd) 

8 

0. 

Pd) 

8 

0. 

PCP(I) 

= 

0. 

PSIA 

8 

.FALSE 

m.ERR 

m 

.0005 

RKT 

8 

.FALSE. 

SI 

8 

.FALSE 

SUBARtt) 

8 

0. 

SUPAR(I) 

8 

0. 

WFL0W 

8 

0. 

XP(I) 

8 

1. 
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$0DK, set in subroutine 0DKINP 


G0NDEL 

= 

l,0E-6 

DEL * 

ss 

.001 

EPS 

es 

0. 

EPSEL 

« 

l.OE-5 

HI 

= 

.01 

HMAX 

s 

.10001 

HMIN 

m 

.005 

IDYSCI 

xs 

0 

IWALL 

« 

1 

JF 

» 

0 

JPFLAG 

» 

0 

JPRNT 


-1 

TEXPLI 

8 

0. 

XMa) 

« 

0. 

XMFTST 

= 

l.OE-3 

$TRANS, set in subroutine TRAN 

ALI 

8 

0. 

IBUG 

8 

0 

MP 

8 

so 

PM GRIT 

S3 
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$TDK (cont'd) 
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